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Der Verdichtungsstoß bei der stationären Umgtrömung 
flacher Profile | 


Von Kl. Oswatitsch in Emmendingen 


. Der Verdichtungsstoß an flachen Profilen verursacht. Entropieanstiege, welche für die Wärbelbildung 
im allgemeinen bedeutungslos sind, die aber dennoch Widerstand bedeuten. Damit hängen eine Anzahl von 
Fragen zusammen, welche nicht immer richtig behandelt wurden und die hier nun geklärt werden. Abschließend 
wird eine Abschätzungsformel für den Ort des Verdichtungsstoßes im lokalen Überschallgebiet gegeben. 


The shock-wave at flat profiles is increasing the entropy, an effect which is generally unimportant with 
respect to the generation of vortices, but is meaning resistance. In connection with that, some questions are 
explained ihat were not always treated. correctly. Finally, an approwimative formula is given for the place 
of the shock-wave in the area of local supersonic flow. 


L’onde de choc sur des profils plats augmente l’entropie. En general cette augmentation de l’entropie 
ne donne pratiquement pas naissance a des tourbillons. Cependant elle repr&sente une certaine resistance. 
Beaucoup d’autres questions sont fonction de ce phenomene. Ces questions qui n’ont pas E£ traitdes d’une 
fagon correcte jusqu’a present sont mises au point. Enfin l’auteur donne une formule d’approximation 
pour L’emplacement du choc de compression dans le domaine suspersonique local. 
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1. Einleitung 

Die Geschwindigkeitsverteilung an einem flachen Profil in ebener Überschallströmung ist 
im wesentlichen bekannt [1]. Die Betrachtungen auf diesem Gebiete können daher keine ent- 
scheidend neuen Erkenntnisse bringen. Es sollen aber einige Unklarheiten geklärt und einige 
im Schrifttum auftretende Irrtümer berichtigt werden. Eingehender wird das Verhalten der 
Strömung in größerem Profilabstand behandelt. Die hier gewonnenen neuen Erkenntnisse sind 
dabei eher für das allgemeine Verständnis als für die technische Praxis von Wert. Für schallnahe 
reine ‚Überschallströmung sind Korrekturen ‚der bisher benutzten Formeln erforderlich. Diese 
entsprechen der größeren Winkelabhängigkeit des Druckes mit Annäherung an die Schallgeschwin- 
digkeit. Eine näherungsweise Berechnung der Verluste im Verdichtungsstoß am Ende eines lokalen 
Überschallgebietes führt zu einer Beziehung zwischen Sprunghöhe im Stoß und Druckverteilung. 


2. Einige bekannte Eigenschaften der stationären Übersehallströmung 
Eine drehungsfreie Gasströmung wird durch das bekannte Gleichungssystem 
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von gasdynamischer Gleichung (1) und Drehungsfreiheit (2) beschrieben. (c ist die örtliche Schall- 
geschwindigkeit; x, y sind die Koordinaten, vw, v die Geschwindigkeitskomponenten.) 

Für schlanke, wenig angestellte Profile läßt sich bis zu hohen Überschallgeschwindigkeiten 
v klein sowohl % wie c gegenüber annehmen 
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Rn Ar noch nicht exakt behandelt. 5. 


 Glauertschen Analogie, in welcher M Hhren den Wert der Mac Hechen Zahl im A instı w > | 


mit®o, als v-Verteilung auf der z-Achse u- 0). DIR Funktion 07 hat das RE :— YM2 M2 in Bei 


Wählt man die Scharen Gleichung (7) und Gleichung (8) als Koordinktenlinien, so vereinfacht 


‘ charakteristischen Grundkurven bei der ebenen drehungsfreien Überschallströmung im Hodo- a | 
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- Gleichung (5) unterscheidet wich ER Ber BR RER, der. ran ‚dt as I, 


gebiet M„ ersetzt wird. re Be 
Für diese linearisierte NH (M= M. >1) It die llgemeine Lösung be- 
- Soll die SERRUNE vor dem Profil ungestört sein, soist y ER $ e Y 
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Sind die Zustände konstant. Es sind dies die Mac h schen Linien, die charakteristischen Grund- DR 


kurven des linearisierten Systems. a 
Neben Gleichung (7) gibt es noch eine zweite Schar von Ma chschen. Linien C5 
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sich ünser Gleichungssystem bedeutend. 
Wählen wir allgemein in unserem exakten Gleichungssystem (A), (2)die Mach schen Linien 
als Koordinatenlinien, so verwandelt sich das Gleichungssystem in [2] 
Ldw_ 1 Fa 
. add Yym-i 
mit w und Vals Betrag und Richtung der Geschwindigkeit und als Machschen Winkel. In 
den beiden Gleichungen (9) entspricht das positive Vorzeichen den Mach schen Linien positiver 
Neigung. i 
Die Bogenlänge kommt in Gleichung (9) nicht vor. Da M oder & nur von w ab- 2 
hängen [3], ist auf einer Machschen Linie 9 nur Funktion von w. Damit liegen die ve | 


graphen unabhängig von den 
Anfangsbedingungen ein für alle- 
mal fest [4]. 


Bild 1. Profil mit stetiger Wandkrümmung in Überschallströmung Bild 2. Stoßpolare und Oharakteristik 
Machlinien; Verdichtungsstoß) 


Alle aus dem Anströmgebiet kommenden Mach schse Linien (Bild 1), welche keinen Ver- 
dichtungsstoß durchkreuzen, entsprechen im Hodographen nach Gleichung (9) einer einzigen 


Die Geschwindigkeit am Profil, Bi, hängt danach RN owned 
_ sprechenden Oberflächenelementes ab, S in Bild 1 fallenden M ac h ac Linien ‚ke 
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u An einer Kreuzungsstelle von Mac h de Linie Rs Verdichtungsstoß ist a ee 
Ste En andie durch die Stoßpolare gegeben [5]. Diese hat die bekannte Eigenschaft, im 
Punkte verschwindender Zustandsänderung mit der entsprechenden Charakteristik Tangente und 


Krümmung gemeinsam zu haben (Bild 2). Die Entwicklung von Stoßpolare und Charakteristik 
AR a kleine nn stimmt also bis einschließlich el Glieder von zweiter. Ordnung überein 
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3. Das Abklingen von Kopt- und Kegel im aiben. ; 
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Wird das Produkt 5 (u, — u) aus (11) und (12) eliminiert, so folgt für schwache Stöße: ch 
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#  eota und cot 5, an der Stelle y entwickelt, ergibt 
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Die Front eines schwachen Verdichtungsstoßes ist die Symmetrale der Machschen Winkel 
“vor und hinter dem Stoß. Die Entwicklung (14) versagt lediglich bei sehr hohen Mach zahlen. n 
Auch das schallnahe Gebiet ist auszuschließen, da dann die Entwicklung (10) schlecht konvergiert, r 
wie später gezeigt wird. 

Es sei nun das Abklingen (der Kopfwelle an einem endlichen Keil untersucht (Bild 3, N 
$.132). Aus den Machschen Linien vor und hinter dem Stoß ergibt sich mit Gleichung (15) 
ein Abklingen. Das Ende des keilförmigen Anstieges sei der Ursprung eines kartesischen Koordi- 
natensystems. Ist die Verdichtungsstoßkurve durch 4 (x) gegeben, so ist 
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Dies in Gleichung (13) eingesetzt ergibt die gewöhnliche Differentialgleichung 
de...1 
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day 2 ( + cot x.) (16) 
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hist der Abstand des Punktes, AN EN - 3 
Stoßes beginnt vom Ursprung des Verdünnungsfächers. 

Die Differenz cot y — cot a, ist ein Maß für die Stärke 
des Stoßes. Nach Gleichung (17) klingt er also mit dr Wur- 
zel seines Abstandes vom Keile ab. Die Verdichtungsstoß- — 
kurve selbst ist 8 durch: ne 
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Bild 3. Be am endlichen Keil 2 
x | Sr BT ra cot 4u — 2 (col &, -estya Y# ... (18). 


(Diese Resultate habe ich im Jahre 1945 in Göttingen vorgetragen. Inzwischen wurde die Auf- 
gabe noch von Lighthill, NPL, London, in einer Arbeit behandelt, die mir augenblicklich 
nicht zugänglich ist.) 

Die Rechnungen setzen voraus, daß der Keil im Punkte = y= 0 sein Diekenmaximum 
hat. Wie der Körper stromabwärts von diesem Punkte fortgesetzt wird, ist gleichgültig. Die 
Machsche Linie der Neigung a, hinter der Kopfwelle sei als „neutrale Machsche Linie“ 
bezeichnet. Ihrer Neigung nähert sich die Verdichtungsstoßneigung mit zunehmendem y. Hin- 
gegen zeigt Gleichung (18) auch,, daß der Abstand des Verdichtungsstoßes von der neutralen 
Mach-Linie mit y über alle Grenzen wächst. Das ganze Strömungsfeld vor der neutralen 
Machschen Linie ist unberührt von dem, was hinter dem Dickenmaximum des Körpers ge- 
schieht. Dies gilt nicht nur für Keile sondern ganz allgemein. Innerhalb der Genauigkeit unserer 
Entwicklungen sind die Ereignisse stromabwärts der neutralen Mach - Linie unberührt von den 
Ereignissen vor dieser, wie wir anschließend sehen werden. 

Zur Berechnung des Verlaufes einer Schwanzwelle am Ende eines sich verjüngenden Keiles 
(Bild 4) seien ganz analoge Größen eingeführt wie bei Bild3. Damit gelten die Gleichung 
(17) und (18) genau so für unser Beispiel der Schwanzwelle. Beide Fälle unterscheiden sich nur 
dadurch, daß bei der Kopfwelle stets y, > a, und daher y >a,, bei der Schwanzwelle hingegen 
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Yo <a, und daher überall y <«, ist. Das Gleichheitszeichen gilt erst für > 0. InBild 3 


und 4 ist die durch den Punkt y= h des Verdichtungsstoßes gehende Stromlinie eingetragen. 
Diese Darstellung des Verdichtungsstoßverlaufes ist also auch für bestimmte gewölbte Profil- 
formen verwendbar, 
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Bild 4. Schwanzwelle hinter dem Bild 5. Abklingen der Kopfwelle an 
endlichen Keil einem vorgegebenen Profil(— — — 
( Stoß; Machlinie; ea ER 


Kopfwelle; 
—— — Stromlinie) ——— Machlinie) 


Hoher Länge. Füı 
amit läßt sich das 
2 ets dur eichung (17) und (18) darstellen, wenn das h nur g 
0% ist, er r Winkel yo ist dabei der Verdichtungsstoßwinkel im Punkte y= h. Entsprechen 
vs, y ai für die 'Schwanzwelle. Das Abklingen von Verdichtungsstößen & an schlanken FriDen ist weit 
a N ‚draußen. also stets der Wurzel aus dem Körperabstand proportional. En ER 
BR:  Konkav. gekrümmte Profilstellen interessieren in diesem ne nike a 
ir ein Verdichtungsstoß 'angefacht. Die Betrachtungen a Ba erst a, des m 
u,“ malen Druckanstieges des Stoßes beginnen. . 
a FE ie Resultate gelten bei Schallnähe nur ungenau. Flache Profile aanon kleine Störunge 
ne Parallelströmung. gewährleisten. Bei schallnaher Anströmung führen aber bereits klein 
' Störungen in die Nähe der Schallgeschwindigkeit, wo nach Bild 2 die Stoßpolare nur ungenau 
durch die Charakteristik dargestellt wird. Je näher M, an 1 ist, desto schlanker muß das Profil 
‚sein, oder desto größer muß die Entfernung sein, für welche Gleichung (17) gilt. Kr 
. Die linearisierte Theorie, Gleichung (6), gibt kein Abklingen i im Unendlichen. Sie stellt als 
dort die tatsächlichen Verhältnisse falsch dar. Dies steht sehr im Gegensatz zur linearisierten 
Theorie im Unterschallbereich, wie sie etwa durch die Prandtlsche Regel gegeben ist. Für 
M, <1 ist eine richtige Darstellung des Unendlichen unbedingt erforderlich. Da die Strömung 
im Unendlichen in allen Richtungen gleichen Charakter hat, kann die Anströmung. nur richtig 
erfaßt werden, wenn das gesamte Abklingen der Störungen richtig dargestellt ist. Für M, >1 Be 
0 sind nur diejenigen Stellen von Bedeutung, von denen aus Machsche Linien auf das Profil auf- 
treffen (in dessen Einflußgebiet das Profil liegt). Was weiter stromabwärts im Endlichen oder 
" Unendlichen geschieht, ist uninteressant. Das Gebiet, welches einen Einfluß auf das Profil aus- 
übt, heiße das Abhängigkeitsgebiet. Je besser dieses i in der linearisierten Theorie BEEN wird, 
desto brauchbarer ist die Lösung. W Na j ii 
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4. Der Stoßwellenwiderstand 
In einer in Göttingen verfaßten Arbeit [7] zeigte ich, daß sich der Widerstand W eines 
stationär fliegenden Körpers durch einen Entropiestrom ausdrücken läßt. Es ist:. 


T, s, esind absolute Temperatur, Entropie der Masseneinheit und Dichte, w, ist die Normalkompo- Pohl 
nente der Geschwindigkeit auf das Flächenelement d/f. Die Integration ist über eine Fläche F zu Some 
erstrecken, welche alle vom Körper verursachten Entropieänderungen umschließt. 


(In der wegen Raummangels sehr kurz gehaltenen Arbeit [7] ist nicht erwähnt, daß die 
E7. Fläche F der Gleichung (1) so groß gewählt sein muß, daß 
r lvdie Reibungsschubspannungen auf ihr verschwinden. Dies 
ist eine übliche Voraussetzung bei allen Impulsmeßverfah- 
ren. Im übrigen läßt sich die Ableitung mit ein wenig mehr 
Aufwand unter Berücksichtigung der Reibungsschubspannun- 
gen durchführen.) 


Schließt man, wie bei allen Betrachtungen dieser Arbeit, 
die Reibungswirkung aus, so gibt es nach (19) dann und nur 
dann einen Widerstand, wenn Entropieänderungen innerhalb 
von F auftreten. Dafür kommen nur Verdichtungsstöße in 
Frage. Also kann es ohne Stoßwellen keinen Widerstand 
geben, andererseits verursachen Stoßwellen stets Widerstand. 
In einer reibungsfreien Strömung kann die Kontrollfläche Bild ee ee 
ohne wesentlichen Fehler unmittelbar entlang beiden Seiten 
der Stoßfront geführt. werden (Bild 6). Die exakte Lösung zeigt allerdings ein Temperatur- 
profil hinter der Schwanzwelle, welches sich durch Wärmeleitung ausgleicht und hiermit zusätz- 
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Y liche Entropieanstiege im Nachlauf ergibt. Diese sind aber gegenüber den Stoßverlusten unbe- a 
Br: deutend, innerhalb der Theorie drehungsfreier Strömung aber überhaupt zu Nr 1% 
# Innerhalb der Näherung drehungsfreier Strömung herrscht nun ‚hinter der Schwanzwelle | 
und vor der Kopfwelle derselbe Zustand 4„, 0x, 7». Das Produkt o w, ist aber auf Deilen SR £ 


des Stoßes gleich. Daher ist / | Sr 
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mit bals Breite des Profiles. Ist As der Entropiesprung im Stoß und wird die Kontrollfläche ganz 
an den Stoß herangerückt, so folgt, aus Gleichung (19) Fer 


W ae 
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wobei das Integral nun über alle Verdichtungsstoßkurven ( zu erstrecken ist. Das zweite Integral 
gilt dabei nur für ein ideales ‘Gas, was hier kaum eine Einschränkung bedeutet (p = der Druck; 
€ €, spezifische Wärmen bei konstantem Druck bzw. konstantem Volumen). 
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Nun ist der Entropieanstieg schwacher Verdichtungsstöße proportional der dritten Potenz 


des Druckanstieges. Man findet mit k = = 
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Der Druckanstieg ist für kleine Ablenkungen wieder proportional der Ablenkung oder auch der 
Größe (Gleichung (17)): 
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In großer Entfernung fällt der Entropieanstieg darnach proportional 
h 3/2 R | 
As = prop (5) ERDE N a | (23). 


Dies in Gleichung (20) eingesetzt und über % integriert gibt, wie erforderlich, einen endlichen Be- 
trag für den Widerstand W/b der Profilbreiteneinheit. Es ist die durch die Verdichtungsstöße in | 
die Luft in Form von Wärme getragene Energie. Keineswegs ist dies ein zum Druckwiderstand | 
hinzutretender Widerstandsanteil. Dieser Stoßwellenwiderstand ist nur die ener- 

getische Ausdrucksform des sich an den Oberflächen in Druckkräften äußernden selben Wider- 

standes. 


Diese Betrachtung lehrt aber, daß die Verluste, welche in der Kopf- und Schwanzwelle 
auftreten, keineswegs unbedeutend sind. Der Entropieanstieg ist in der Flächeneinheit wohl 
gering, erstreckt sich aber über große Flächen. 


Je flacher das Profil, desto größere Werte nimmt % in Gleichung (17) an, desto langsamer 
nähert sich also die Richtung des Verdichtungsstoßes dem Winkel &,. Bei flachen Profilen ist 
dabei ein großes Gebiet für den Stoßwellenwiderstand zu berücksichtigen. Dielinearisierte Theorie(6) 
schließlich entspricht einem verschwindend kleinen Entropieanstieg im Stoß auf unendlich 
langer Strecke. So ist es zu erklären, daß scheinbar im Gegensatz zum Integralsatz (19) oder (20) 
mit der inearisierten Theorie Druckwiderstände errechnet werden. 


Dies zeigt, daß es sinnlos ist zu fragen, ob der Entropieanstieg vernachlässigbar ist oder 
nicht. Dies hängt nicht einmal davon ab, was man berechnen will, sondern wie man die Rechnung 
durchführt. Man kann den Entropieanstieg vernachlässigen, wenn man den Widerstand flacher 
Profile aus der Druckverteilung am Körper berechnet, man kann es natürlich nicht tun, wenn man 
den Weg über den Stoßwellenwiderstand geht. Unterschallströmungen mit lokalen Überschall- 
gebieten müssen auch Verdichtungsstöße aufweisen, wenn Druckwiderstände auftreten. Nähert 
man den Strömungsverlauf nicht grade durch unendlich hohe Überschallgebiete, so muß ein Wider- 
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=; ist die in üblicher "Weise: definierte Stromfunktion. ' Die Entropie der Masseneinheit, ist Er 
‚ Stromlinien konstant, ‚hängt also lediglich von Y ab [8]. HR | 2 ee F 


' Für die Größe der auftretenden Wirbel ist nach ( (24) also die Entropieänderung quer zur N 
SER ruhe maßgebend. Hier kommt es auf das Verhältnis von Entropieanstieg im Stoß zur Länge 
des Verdichtungsstoßes an. Da das Produkt von Entropieänderung und Stoßfrontlänge wegen dr 
—_  Endlichkeit des Widerstandes. sicher endlich ist, ist bei flachen Profilen, d.h. schwachen aber 
en Stoßwellen, das Verhältnis von Entropieanstieg zur Stoßlänge klein von der zweiten ‚Ord-" re 
nung. Dies ist der physikalische Grund, weshalb nicht nur in der linearisierten Theorie, sondern ; 
auch noch in einer nächsten Näherung drehungsfrei gerechnet werden darf. N Ri. 
©... Die Integration von Gleichung (2 ) führt zum bekannten Integralsatz 
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wyist die Tangentialkomponente der Geschwindigkeit w an das Längenelement dlder geschlossenen f 
Kurve (. P 
Br ‚Als Kurve c sei die aus dem Unendlichen kommende, am Profil unmittelbar vorbeiführende . 
Stromlinie gewählt und durch zwei Gerade senkrecht zur Strömung und eing Stromlinie in Fo re 

- Entfernung zu einer geschlossenen Kurve ergänzt (Bild 7). Da die götromimie in rer Bao 
keiten, wie die Potentialtheorie lehrt, mit dem Quadrat der Entfernung abnehmen, die Kurven- 
länge von 1 aber nur mit der ersten Potenz’der Entfernung wächst, läßt sich das Integral Gleichung 
(25) auf ein Integral über die Körperstromlinie vereinfachen: 
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Werden die Geschwindigkeiten über die Bogenlänge der Stromlinie aufgetragen, so sind die von 
den Über- und Untergeschwindigkeiten gebildeten Flächen gleich groß (Bild 7). Dies gilt auch 


für kompressible Strömungen, selbst wenn lokale Überschallgebiete eingebettet sind, solange nur > DS; 
! drehungsfrei gerechnet werden darf. Auch dann klingen nämlich die Störungen in großer Ent- 3 
z fernung mit dem Quadrat ab, wobei lediglich der Faktor 1/) 1— 12 auftritt. Dies stört aber nicht _ e 
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Bild 7. Die Zirkulation an Profilen Bild 8. Abklingen der Störung in der Überschall- 
in Unterschallströmung strömung in großer Entfernung 


In der linearisierten Überschalltheorie hat es nun keinen Zweck, eine in das. Unendliche 


reichende Kurve zu ziehen, da im Unendlichen das gleiche passiert wie am Profil selbst. Doch 
folgt allein schon aus der Gleichung (6) wegen [ "de=0 
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keiten vor dem Dickenmaximum und den Überg 


 Unendlichen ab. Das Verhalten in großen Entfernungen sei daher näher untersucht. 
© Während sich ein Körper in Unterschallströmung in großer Entfernung stets wie ein Dipol = 
" verhält, kann das Verhalten eines Körpers in Überschallströmung näherungsweise durch die 
‚Strömung an einem Rhombus dargestellt werden. Weit ab erscheint die Expansion an einem 
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rebildeten Flächen sind einander ech [00 
en ara N ärehungsfreien Überschalifeld mit Tilfe von Gleichung (9) 
die genauere Geschwindigkeitsverteilung rechnen, so würde man sehen, daß Gleichung (26) 
nicht mehr erfüllt ist. Das erscheint zunächst verwunderlich, denn vor der Kopfwelle und hint 
der Schwanzwelle herrscht ungestörte Parallelströmung (u„) und die Störung selbst klingt 
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Körper stets wie die Prandtl-Meyer- Strömung an einer Ecke (Bild 8). 
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ee | 
u: “ 


Die Abszisse der neutralen Mach -Linie im Abstand 9, sei z,. Der Machsche Winkel | 


im Abstand y, ist dann gegeben durch Se 


% In | j ; 
cota= ul ur, Erz 2 RAAB 


Hn n 
Die Störung der Geschwindigkeit in weiter Entfernung kann in erster Näherung ersetzt werden 
durch: 


u = () era ta) +... age EN a 
2 M*:coto, ([C—%, 
ER ER M: =( Yn )+-- 


M* ist die durch die kritische Schallgeschwindigkeit c* dividierte Geschwindigkeit w. Gleichung (29) 
kann nach bekannten Formeln [3] leicht gefunden werden. 

Die Grenzen der Störung auf der Geraden y— 4, sind durch Gleichung (18) für Kopf- und 
Schwanzwelle gegeben. Es sei die Neigung der Kopfwelle im Punkte y= h: yo ,;, jene der Schwanz- 
welle: Yg,, dann sind deren Abszissen im Abstand %,„ gegeben durch Gleichung (18) 


I = %n — 2 (Cot &o — ot Yoz) Y nf Nr (30). 
= %y +2 (cotyos — cot aa) Y Yahı 


[4 


Das Integral-über die Störungen im Abstand y, ist damit: 


% u 
ZN c—% 
fre—un as ung gr otan | de > 
k+1 M& 
677 a 17 In we) (31). 


2. Mk: 
= 2 Us k-ıMe cot &% h [(cot X — cot Yo)? — (Cot yy, — cot &,)°] 


Durch A, yoz und yo, ist die Profilform gegeben. Der Ausdruck in den eckigen Klammern enthält 
die beiden unabhängigen Winkel yo; und yo, und ist daher im allgemeinen #0. Wesentlich ist, 
daß die Integration einen von y, unabhängigen Ausdruck ergibt. Das Integral verschwindet also 
keineswegs in großer Entfernung. Damit ist gezeigt, daß Gleichung (26) bei Überschallströmung _ 
exakt nicht erfüllt ist. Das Integral (26) kann positive oder negative Werte annehmen und 
kann auch verschwinden. 

Wenn nur vorausgesetzt wird, daß die Strömung isoenergetisch ist, daß also über die Strom- 
linien hinweg keine Energie fließt, so bleibt die gasdynamische Gleichung (1) unverändert. Die 
Entropieänderungen in den Stoßwellen verursachen lediglich eine Änderung der Gleichung (2) 
für die Drehungsfreiheit, welche-nun durch Gleichung (24) zu ersetzen ist. Die Machschen 
Linien als Koordinaten eingeführt, ergeben nun die Gleichung [9] 

er dw  cota& ds 

FdP-+cota 2 See ae a ee A Be (32). 
Das obere Vorzeichen kommt den steigenden, das untere den fallenden Machschen Linien zu. 
Zur entsprechenden Gleichung (9) ist hier ein Glied hinzugekommen, welches die Entropieände- 
rungen berücksichtigt. Wenn sich diese quer zur Strömungsrichtung ändert, so liegen die 


to ‚sind N Bon bilden re 
e Ma welche aus, dem. Nachlauf fallend. ‚ode 
r Kpfwe in das (ine Berügen wird, ‚ist 49 = =0. N 


Br % 
Br x ER; pr = Ks 
 Gichung ee 1äßt Re exakt RR ‚Die größten a ofen: an len 
Stellen höchsten Entropieanstieges,. ‚also im allgemeinen unmittelbar an der Profilstromlinie 


Nach Gleichung (22) steigt die Entropie mit der dritten Potenz des Drucksprunges im V 
Fa ongsstnE an. ‚Der Entropieanstieg nach an ee a oder nach dem . i 


flachen Profilen abe Heiden sich also die Geschwindigkeiten im Nachlauf - — wie, erwartet - 

kaum von der Anströmgeschwindigkeit. "Wenn allerdings die Ablenkungen’und damit die Entropie- 
- . „anstiege in der Schwanzwelle erheblich sind, so können starke Nachlaufeffekte auftreten, obwohl. 
x." zwischen ‚Kopf- und Schwanzwelle drehungsfrei gerechnet werden darf. Bei der ne 
| ee rag am Keouß sind starke a a sogar die ee SENT 


AB; Die Berücksichtigung der Entropieänderungen hei der Druckverteilungsberechnung s 
‘Um das erste durch die Entropieunterschiede bedingte Glied zu berechnen, genügt es, in. 

“ Gleichung (32) den Koeffizienten des letzten Gliedes als konstant gleich der entsprechenden Größe ' 

- © “in der Anströmung zu setzen. Gleichung (32) Ian sich dann Nihgerionen. Für. die fallende Mac h - 

: 

3 

| 


"Linie gilt: 
2 


Er os (S3—- 54) 


b-M2 G—% 


h %,—9,= | cot a + (34). | 
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Damit ist der Entropieeinfluß auf der Machschen Linie durch die Entropie im Anfangs- und 

Endpunkt ausgedrückt. Die Geschwindigkeit nach dem Winkel entwickelt BR eine Reihe von 

der Form: | 
wu 1 89,8 


4 —d)) +0, (da — dB) +65 OBEN el ei en 85): 


£ Um den Druck nun längseiner M ach schen Linie zu entwickeln, kann nicht mehr der isentropische 
. Zusammenhang von Geschwindigkeit und Druck zugrunde gelegt werden. Es ist die Energie- 
gleichung zu nehmen, welche in einer in einer früheren Arbeit 17] gezeigten Weise geschrieben 

werden kann: BR 


\ i j 
wdw + dp +Tds=0 ar BERGER. 1136); 


: Gleichung (36) können alle geforderten Ableitungen für die Funktion p (w, s) entnommen werden. 
Man erhält folgende Entwicklung: 


1 /®p 1 (8®p\ - u 
Pa Pı = — 01%ı (wg — W)) no (=) (w, — wı)? aer® (5) (wW-wP-T, 01a) +: +: 
E = (37). 


Hier unterscheiden sich die ersten drei Glieder nicht in der Entwicklung bei isentropischer Zustands- 
änderung. Wir haben nun Gleichung (35) in (37) einzusetzen, um die Entwicklung der Druck- 
änderungen nach den Winkeländerungen und nach der Entropieänderung zu erhalten. Bei Ver- 
nachlässigung aller höheren Glieder geht die Entropie durch Gleichung (35) nur in das Glied erster 
Ordnung in (37) ein und hebt sich mit dem letzten Gliede gerade wieder auf, Die Entwicklung 
des Druckes nach dem Strömungswinkel auf der Charakteristik ist danach bis zu Gliedern 


a A a a lat a = 25 


138 \ Oswatitsch, Der Verdichtungsstoß bei der 


von einschließlich dritter Ordnung unabhängig vom Entropieverla . Es sei ER 

m  Mernd +, tar een (8). a de 
1 T £ Nie 

| 2 N Kaas 

Für die Entwicklung der Druckänderung mit der Winkeländerung im Verdichtungsstoß schreibt RX 

man. analog [10]: % 

| 1 Belae GI) + HD) d—d)t.:. - - 


2 
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Es sei nun mit dem Index 1 die Anströmung gekennzeichnet. Zur Berechnung des Druckes p, 

"an einer’ Stelle mit dem Winkel 9, am Profil ist die Ma ch’sche Linie stromaufwärts durch die 
Kopfwelle hindurch in das Anströmgebiet zu verfolgen. Die gesamte Winkeländerung setzt sich. 
nun aus der Winkeländerung im Verdichtungsstoß an der Kreuzungsstelle mit der genannten 
Machlinie und aus der Winkeländerung längs dieser zusammen. Diese beiden Winkeländerungen 

. aber haben nach Gleichung (38) und (39) verschiedene Wirkungen, weshalb es allgemein nicht 
möglich ist, p, bis zu Gliedern dritter Ordnung von 9,—, zu berechnen ohne die Ablenkung von 
Verdichtungsstößen am Ausgangspunkt der Charakteristik zu kennen (Bild 1). Danach ist also 
zur Berechnung der Druckverteilung mit der geforderten Genauigkeit eine Konstruktion erforder- . 
lich, die im wesentlichen auf eine Charakteristikenmethode hinausläuft. : 

:In den bisher angegebenen Formeln wurde einfach die Ablenkung der Strömung an der 
Profilspitze genommen. Die Zulässigkeit dieser Vereinfachung hängt davon ab, wie weit sich die 
Winkeländerungen am Profil hinter der Kopfwelle auf die Strömung am Profil weiter stromab 
auswirken können. Maßgebend dafür ist, wo die von der Körperspitze ausgehenden und an der | 
Kopfwelle reflektierten Mac h schen Linien wieder auf das Profil auftreffen. Diese Betrachtungen T 
lassen aber den Aufwand, welche diese dritten Näherungen verursachen, kaum gerechtfertigt er- M 
scheinen. ‘Bei geringeren Überschallgeschwindigkeiten oder starken Profilkrümmungen in der } 
Nähe der Profilspitze sind die genannten Vereinfachungen jedenfalls falsch. 


7. Die Geschwindigkeitsverteilung bei schallnaher Überschallströmung 
Für die Betrachtungen dieses Abschnittes sei stets drehungsfreie Strömung vorausgesetzt. 
Bereits bei schallnaher Strömung sind die Entropieanstiege im Stoß gering, die Stoßfronten aber 
lang, so daß gerade hier die Annahme der Drehungsfreiheit dem wirklichen Verhältnis weitgehend 
gerecht wird. Die Entwicklung des Strömungswinkels nach der Geschwindigkeit auf-der nun | 
im Hodographen festliegenden Charakteristik im Punkte (9 = 0, w= w,) kann in folgende Form | 
gebracht werden: 
Mh N ne 224 n 
RLENEPR N EA z{i aaa mi) sel. 3 M en 
(M2— 1? Mm? —1 ILL TOR mı) r 6) Mm-—ı SET: 


Es ist dies die Umkehrung von Gleichung (35) bei Isentropie. f(M,) ist eine Funktion der Mach- 
schen Zahl, welche bei Schallgeschwindigkeit den Wert: 

| k-+1\2 

)=-— ( >) 


ad 


? 


da GR N A oe 


annimmt, 
Für Schallgeschwindigkeit (w= c*) wird nun der Ausdruck 


Le! u 

a ES PD EHE 
M-—ı M,-—-ı)ıM, +1) 
Die Entwicklung (40) konvergiert damit nahe der Schallgeschwindigkeit ziemlich schlecht. Daß 
Stoßpolare und Charakteristik in beiden ersten Gliedern übereinstimmen, ergibt noch keine sehr 
gute Übereinstimmung beider Kurven bei Schallgeschwindigkeit. Die Charakteristik verzweigt 
sich eben an der Schallkurve, und diese Tatsache ist durch eine Potenzreihe der Form (40) nicht 


0,5. 


für Fe Stremung wohl hl am besten e: 
wicklung welche wie e folgt beginnt: | 


Zum ER Vergleich hie Ar Stoßpolaren s sei die Entwicklung 5 in ı den Geschwindigkeitskompo- 
 nenten Bestehen. Dabei BELLE, DE SEEN \ 
- \ a REN für vu EN. ER ER (43). 
Nach einiger ah u man: | ar BE I N ie 


/ 


\n—cr 


ken man er die a si für kleine Abweichung von der Schatlgeschwindigkeit, 2 her. 
so findet man nach einiger Rechnung. für die Glieder erster Ordnung Nr 


% Dean Y”y or! en uch: ar = k N 
2-5 (a lese er a) 


DE RR OR: 0). 
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2 Damit sind Ar Ansdriieke unter der Wurzel alle in derselben Erößenertunue, ‚Bei kleinen Über Br. e: 
“ schallgeschwindigkeiten von %, unterscheiden sich also N v von ER und Charakteristik 
r fur u = c*: 
E07 Le /2 Gas, 3/2 0 oe 
EEE 3 &- ) k+1 > Al 
Er CE ARE R 1— j* 1 — 0,057 | Br rs 1) = 2 Naalie AZIE 
E | Während die Charakteristik (44) senkrecht in die Schatlinie (in unserer Näherung u = c*) mündet, 
7 wird der größte Wert von v/c* der Stoßpolaren bei Unterschallgeschwindigkeit 


3 | Ya 2 Er EN 1 
= a N ee eng (a1) bei ee EN Eh (48) 


- erreicht. Dies N in unserer Naheraüp der maximal zulässigen Ablenkung bei Anliegender Br, r 
Kopfwelle. KERN 
In, der Praxis werden im allgemeinen Druckverteilungen an Profilen von etwa 10% Dickes Yo 
gerechnet. Bei diesen wird die Schallgeschwindigkeit an der Profilspitze etwa bei einer Anströmung ni 
mit 1,3facher Schallgeschwindigkeit erreicht. Ein kleiner Fehler, welcher aus der Verschiedenheit BR 
- von Charakteristik und Stoßpolare resultiert, muß wohl hingenommen werden, wenn man Ge- Be: 
schwindigkeit oder Druck als Funktion der Stromlinienneigung allein darstellen will. Der Fehler, Bu 
welcher aber in der schlechten Konvergenz der Entwicklung (40) begründet ist, ist zu umgehen, 
indem man bei schallnaher Strömung mit einem Glied entsprechend Gleichung (42) beginnt. & 
Allenfalls wäre es am besten, mit einer Kurve zu arbeiten, welche zwischen Stoßpolare und Cha- 
rakteristik liegt, um auch dem Kopfwelleneinfluß etwas gerecht zu werden. 
Wegen der höheren Winkelempfindlichkeit der schallnahen Überschallströmung sind hier 
die Druckanstiege höher als bei großen Überschallgeschwindigkeiten. Gegenüber der linearisierten 
Theorie und der Theorie zweiter Näherung [6] sind also höhere Widerstände und größere Auftriebe 
-zu erwarten. (Der letztere Effekt hat ein Gegenstück bei der schallnahen reinen Unterschall- 
strömung, wo. die genaueren Rechnungen auch eine Auftriebserhöhung der Prandtlschen 
Näherung gegenüber ergeben.) Das Profil geringsten Widerstandes hat mehr als 50% Dicken- 
rücklage. 


na 9 Al ar ain u. & ha 
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8. Der Verdichtungsstoß im lokalen Überschallgebiet 
Da sich der Widerstand sowohl durch die Druckkräfte auf das Profil als auch energetisch 
durch die Verluste im Verdichtungsstoß darstellen läßt, kann eine Näherungsformel für den Ort 
des Stoßes gefunden werden. A. Betz hat vor einigen Jahren in einer kurzen, mir augenblicklich 
nicht zugänglichen Arbeit den Ort des Stoßes dadurch festzulegen versucht, daß durch ihn kein 
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| Schub verursacht wird. Während bei A. Betz EN N 


...0......Ortdes Stoßes geschlossen wird, soll hier dt Widerstand st durch Or und öhe des Druc KR 
RR sprunges am Profil ausgedrückt werden. 1 
ER a IS DIET Widerstand W eines Profilstückes der Breite bist TE atee ! 
Br “7 .r — FRI Pe le & .. (EU “ , 
Bee eh: ' - W b di pdy ne RE 6 N % ie, 4 
2 3 Es ist dabei über die Profillinie zu integrieren. Der Einfachheit halber seien die Betrachtungen 
0 auf die obere Hälfte eines symmetrischen nicht angestellten Profiles beschränkt. | 
ER “ nung des -Stoßwellenwiderstandes sei im wesentlichen auf das Hauptglied beschränkt, da FR Be. 
n 2 Genauigkeitlin den folgenden Rechnungen notwendigerweise nur gering ist. Kuppeerben ” 2; 
Da - den Rechnungen in Abschnitt 4 setzt man in Gleichung (19) ‚ ER | 
P, em af= e* w* bdy— Qu un bdy, a 
BEN . woraussich wieder Gleichung (20) ergibt. In diese Gleichung (22) eingesetzt führt zu Gleichung (50): "s 
. k+1 ea FI 
Wen | le Ay. en en (MD). N 


Im Nenner des Integranden schreibt man werkaniBig den kritischen Druck p*, der der mittlere 
Druck im Stoß ist. Da der Druck dort in der dritten Potenz auftritt, seien hier keine unnotwendigen 
j Vernachlässigungen eingeführt. 
x, Zur Berechnung des Integrals der Gleichung (50) ist vor allem die Kenntnis des Druckes in 
Profilnähe erforderlich. Weiter draußen nehmen die Verluste stark ab. Der Druckabfall am 
Profil quer zur Strömungsrichtung ist durch die Zentrifugalkräfte mit R(z) als Wandkrümmungs- 
- radius gegeben. Es ist: £ 


r 


—_ nn m 


nn 


. 1). 


Begnügt man sich mit der groben Näherung eines linearen Druckabfalles bis zum kritischen 
Druck, wobei in Gleichung (51) am besten M=1 und p= p* zu setzen sein wird, so ist: 


ry= P(O) +kp* 5, 

wobei y vom Profil an gezählt wird. 
Schon wegen dieser rohen linearen Näherung, aber auch wegen Verwendung von Gleichung 

(22) für den Entropieanstieg im Stoß, welche höchstens bis M = 1,2 eine brauchbare Näherung 

darstellt, ist es erforderlich, sich auf kleine Überschallgeschwindigkeiten zu beschränken. Für 


solche ist der kritische Druck p* das arithmetrische Mittel der Drucke vor und hinter d&m senk- m 
rechten Verdichtungsstoß und folglich: 4 


‚Ap_ „Pro p(0) Far | 
Fanart al 1]. 
wenn p.(0) der Druck vor dem Stoß ist. Dies in Gleichung (50) eingesetzt und bis zum Stoßende 


?(y) = p* integriert, ergibt: 
iR 1 p 4 
=.d4 Be 
h eur R(' ) 
Dies ist nun dem Wert von Gleichung (49) BR woraus folgt: 


FR Az We ERDE En 
Fe Au ırlı FRE AN TE (52). 


- Der Index s deutet an, daß die’ en Werte an der Stoßstelle zu nehmen sind, wobei 


es in Gleichung (52) gleichgültig geworden ist, ob für p,.der Wert vor oder hinter dem Stoß ge- 
wählt wird. Da für jedes endliche Profil 


ist, kann im Integral (52) eine beliebige konstante Größe zum Integranden addiert werden. Besser 1 
ist es noch, einen Druckverlauf zu subtrahieren, welcher den Widerstand Null ergibt, also etwa 

die mit der Prandtlschen Regel gerechnete Druckverteilung. Dann braucht nämlich nur 

über ein kleines Gebiet von Druckunterschieden integriert zu werden. 
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ing die en Macher Profile zu betrachten. Die San ae am Profil Wirkk e in ic 
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' Die Folge ist, daß das en ‚über die REN in a a 
einen konstanten Wert annimmt, Die über der Bogenlänge der Profilstromlinie AUISCHABehe 
Unter- und Überschallgeschwindigkeit bilden daher nicht mehr wie in der gr oder noch _ 
in der linearisierten Überschallströmung gleich große Flächen. 

Die Entropieanstiege im Stoß sind zwar gering, doch über die langen, URAR. abklingenden 
- Stoßfronten integriert, I. sie endliche Verluste, welche zum bekannten Widerstand in Über- 
schallströmung führen. Für die. Wirbelbildung sind die Entropieanstiege aber meist unbedeutend, 


. weil es dabei nicht wie beim Widerstand auf das Produkt aus geringem Entropieanstieg auf langer 


Strecke sondern auf Änderung eines geringen Entropieanstieges auf großem Abstand ankommt.‘ 


Das Zusammenfallen von Charakteristik und Stoßpolare am Hodographen in großer Um- AR: 


 gebung des Anströmungszustandes lassen die Bestimmung von Druck und Geschwindigkeit am 
- Profil lediglich als Funktion der entsprechenden Stromlinienwinkel innerhalb der Näherung einer 
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zweiten Ordnung zu. Eine Näherung der dritten Ordnung dieser bekannten Formeln ist aber im 


allgemeinen nicht möglich, da zu diesem Zweck der Kopfwellenverlauf innerhalb von einem be- 
stimmten Abstand vom Profil bekannt sein muß. Die ke des Entropieanstieges am Profil 
selbst genügt allein nicht. 

Die Konvergenz der Entwicklungen von Charakteristik und Stoßpolare im Hodoprapiek 
im aaa nisfand wird mit Annäherung an die Schallgeschwindigkeit schlechter. Auch bei 


den schlanksten in Frage kommenden Profilen weicht die Stoßpolare von der Charakteristik in 


der Umgebung der Schallgeschwindigkeit etwas ab. Die mangelnde Konvergenz der Entwicklungen 


- kann hingegen umgangen werden, indem man-auf gebrochene Potenzen übergeht. Dadurch kann 


die größere Winkelempfindlichkeit bei Schallnähe dargestellt werden. Praktisch ergeben sich 
daraus höhere Widerstände und’ Auftriebe und stärkere Dickenrücklagen für Profile kleinsten 
Widerstandes. / 

Für die Druckverteilung an Profilen mit lokalen Überschallgebieten läßt sich eine Integral- 
bedingung ableiten, indem man den Widerstand sowohl aus der Druckverteilung als auch aus den 
Verlusten im Verdichtungsstoß berechnet. Diese muß neben der — wohl schon bekannten — aus 


‚der Drehungsfreiheit folgenden Integralbedingung, daß die an den Über- und Unterschallgeschwin- 


digkeiten gebildeten Flächen einander gleich sein müssen, erfüllt sein. 
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Kontinuumsmechanik angeknüpft werden. Sind ß 
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7 Formänderungsgeometrie der FR: 


Wie bei Herleitung der Gleichgewichtsbedingungen der Schale, so ko ädeh bei 
der geometrischen Bedingungen des ‚Formänderungsvorganges an die 


die kartesischen Koordinaten des deformierten Körpers, so stellen die Größen VF= Y, die 
kartesischen Komponenten des Verschiebungsvektors dar, der als kleine Größe 


vorausgesetzt 
wird. Das krummlinige Koordinatensystem «* möge an die Körperpunkte gebunden BER ar Ze | 


die Peformefionaiewegung mitmachen. An Stelle von Gl. (1) folgt dann: 


rc +Pell. en de 


Mit le auf Gl.(4) ergeben sich als kovariante Komponenten des metrischen Tensors des 
deformierten Körpers unter Beschränkung auf Glieder erster Ordnung hinichtlich vr: 


Kun = Wale a Kalle N 


Die halbe Differenz der kovarianten Komponenten des metrischen Tensors vor und nach der 


Baosmation liefert den Dehnungstensor: 


1 
an 5 (u u) N ORTE rare RE (103). 


Aus Gl]. (102) folgt mithin (vgl. hierüber, sowie über weitere Beziehungen des deformierten 
Körpers, die in Fußnote 8 zitierte Arbeit des Verfassers): 


1 
„= Wiltrld. 0. ern nn. 104). 
Die kontravarianten Komponenten des Dehnungstensors sind entsprechend: 
1 - 
een: U u EN le (105), 
bzw. nach Ausrechnung mittels der ee des zweiten Abschnittes: 
d oVe Agiu 
ER De an 
dan — le == hir 5 v”) a ee . (106). 


Der so definierte Dehnungstensor steht auf Grund des Hookeschen Elastizitätsgesetzes in 
linearem Zusammenhang mit dem Spannungstensor. In der hier verwendeten Schreibweise 
gilt für homogenisotrope Körper: 


Su 2 oa + nn grau d) x Er we ae DESE na (107). 


Hierin ist G der Schubelastizitätsmodul, 1/m die Poissonsche Konstante und 


dm Pr dh VE A Ne (108) 
die Summe der Dehnungen. 


Bei der Schale sind erfahrungsgemäß die Querkraftschubspannungen und die senkrecht 
zur Schalenmittelfläche wirkenden Normalspannungen sehr klein gegenüber den Längsspannun- 
gen und den Biegespannungen der Schale. Bei Anwendung der Formänderungsgleichungen (107) 
auf das Schalenproblem können daher die Spannungskomponenten mit dem Index 3 als ver- 
nachlässigbar klein angesehen werden, so daß aus Gl. (107) folgen: 


BB N I ee © (109) 
und 
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ee 


Kneipe > 
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Wegen GI. (108) und (105) geht GI. (110) über in. 


Fi er‘ a 1 
an le 1 ua 7 Volle Er | (111), 
so daß 
1 1 
er BE EA er ee Beute (112) 


folgt und Gl. (108) für die Flächenkomponenten des ee der Schale die Beziehung 


20 (dr + — —#4)) A a IT 


liefert. Aufgabe der weiteren Rechnung ist die Zurückführung der Komponenten des Dehnungs- 
tensors auf die Verschiebungskomponenten ie und Ir = V; der Schalenmitteltläche und ihre 


Flächenableitungen. Zu diesem Zweck seien REchSE die aA en der Verschiebungs- 
komponenten näher untersucht. Aus Gl. (109) folgt mit Bezug auf die Definitionsgleichung (105): 


Vale rel N SL 
oder wegen Gl. (106): 
. oV* 
=, UBS 
SE a are (115) 


Während die bisherigen Beziehungen dieses Abschnittes noch für alle Punkte der Schale, also 
auch außerhalb der Schalenmittelfläche Gültigkeit haben, bleiben die weiteren Rechnungen 
auf die Schalenmittelfläche beschränkt, worauf der Index 0 hinweisen möge. Für die erste 
Raumableitung der Verschiebungskomponenten der Schalenmittelfläche ergeben sich die maß- 
geblichen Beziehungen aus Gl]. (63), (111) und (114). Werden zur Abkürzung für die Normal- 
komponente der Verschiebung der Schalenmittelfläche und ihre Raumableitungen nach x? die 
Bezeichnungen 


ra==.W, | =D: an SET NER ER RE (116) 
) 0 ||s 
eingeführt, so folgen: 
V»®— V% 
DE 0 
veB—= — lez.dr Vf, 
0 () 
EN N EREN WER 1. (117) 
v2 = W|° +b5 vs 
0 0 
1 
A ey REN N =D 
% | Mel "X & D KR) 
Für die Flächendehnung ergibt sich demnach: 
Mae VA BEN ea es a (118). 
0 & 0 0 .|& 


Die zweite Raumableitung wird unter Verwendung dieser Beziehungen an Hand der Gleichungen 
(63) bis (66) errechnet, wobei sich ergibt: 


Veljß EN ber Fr b% Sg +bP Vay — ca W—bePD, 

o ls er, Bo a) 

rar END; 

v* Ta “— — D|®, 

0 lan o |ls (119) 
y3 “= Wed + b2je Yr + ba Vriß dB Vria — baBD, g n 
0 rg 9 ro 

y33%— D|e, 

Se ’ 

ri eat a) Val 

0 m—1 o.||sa 


> RR, .. WE 
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ee. Nonkehflassenisich VI. BE N: 
BER RANUIRDIeleungeR yes ri ® Mit Bezug auf or 


Hierbei ist | BR ee > 


ein schiefsymmetrischer Tensor, welcher der örtlichen Dh um die Flächennormale an 
SEE bckape gilt | ö re | 


FE | 
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Far die erste Raumableitung des Dehnungstensors nach =° ergibt sich mit Bezug nit am R 


=, N beaD— cr W Fr a AO 


Fyü 


r 


P- 3 ern )- ER E Es 


ah Mes, PuD—aM.... 0... ia 


Für die zweite Raumableitung des Dehnungstensors nach =? folgt aus Gl. (120): [ Era; 
RE EN BEE HER | a See B ya _ 
En ee ee ai 
bl Hu a av +es Der: +4 Beer W—2c#D, . (125). 
de = —Di% —bEF — 20: D— 2bop wi? — 25x Bi], el 
0° 1133 


Nunmehr sind alle in den Spannungs-Dehnungs-Gleichungen auftretenden Größen auf die 
Verschiebungen der Schalenmittelfläche und ihre Flächenableitungen zurückgeführt. Werden 
noch in Analogie zu Gl. (92) für den Dehnungstensor und seine Raumableitungen nach 2° Ab- 
kürzungen eingeführt, indem 


en ran me 


an une ne a 
"di er 8 Er 8 
0 +5 0 or 3 & u 


dab 
0 


= kaß, kx =k | 

# 

gesetzt wird, so gilt 
dl = -- Dief — ber F +53 KPr aa erg + co dar 


0.1133 =7 | 
dl = — Die —bEF +26,pKaP . (127). 
0 ||33 


Zufolge Gl. (92) und (113) gilt weiter: 


Dr =2an(en 4 { 


TE re ar een Sn 


. . 
— ——n nn nn te Seen rn 


1 1 
Me — a ae +; PR) re Bünde; (129). 


Damit sind alle Spannungsgrößen auf die drei Verschiebungskomponenten der Schale zurück- 
geführt, für deren Bestimmung die bisher abgeleiteten Beziehungen, zusammen mit den im 
nächsten Abschnitt behandelten Randbedingungen gerade ausreichen. 


8. Die Randbedingungen der Schale 


Zur Untersuchung der Randbedingungen werden zweckmäßig für Tangente und Normale 
des Randes besondere Bezeichnungen eingeführt. Ist ds Linienelement längs des Randes, 


r EN Ba EEE N En A R% EEE ET 
"die Ben ae er a längs‘ des Ran es. "Der Binheitsektor 
_ Randnoranale ae a EEE auf a. a en * 


Apfıı ; = x 
PR, elle VE > AR 


H 


_ Wofern Sn die nee auf die a Er 'Schalenrandes Deren etsia} 
 rraltanbe) sind die Verschiebungen ” w und die Drehung & des Schalenrandes um 


> Randtangente. vorgeschrieben, owohel sich 17) entsprechend gr ‚17 aus A NINE 


: D—- g% € Re ers —. m oo > RR Ta “ be m. IHRER 3 (132) 
een läßt. ; . ; WAR RN A ER N E 
EN Ur ie Beziehen sich andererseits die ee auf den LE der Schalen t 


(zweite Randwertaufgabe), so sind die resultierenden Kräfte und Momente en % 
welche den längs des Randes wirkenden Spannungen entsprechen. Zu ihrer Berechnung ist die 


Bezugnahme auf die physikalischen Spannungskomponenten 8?# V9uu|9* erforderlich 4). Die 
ZURSROFIEER Komponenten in Richtung =* werden durch Multiplikation mit dem Richtung : 


kosinus V9«„ erhalten, während sich die Kraft, welche an dem von den Differentialen d.x” und. 


Br rar KEN gebildeten Flächenelement wirkt, durch Multiplikation mit Vag* derdze 
errechnet und mithin den Betrag 8?* /gc* da’ dx annimmt. Für die Schale gilt für Yg ent- B3 | 
En: el (Gl. Io) (42) und (43) bei Entwicklung bis einschließlich «2 der Ausdruck: 

: W= ja +2wH HR) RL Pe Rn . (133). 

Z Mit end des in Gl. (35) definierten schiefsymmetrischen Flächentensors kann demnach 

= die auf das Längenelement des Randes bezogene Kraft in der invarianten Form Fi 

2 a | 
Er Ber Przesefsene FIoH ar) a NN 


er 
geschrieben werden. Die Integration liefert mit a auf die Entwicklungen (45) und 


(70), sowie die Gl. (85), (92) und en unter Memeehlissiaune eines mit dem Faktor A?/12 be- 
hafteten Restes: 


Ä je Pk = e,grP [(Lar + Mer + MON EHONN. 2.2... (135). 
Mithin ergibt sich für die tangential zum Rand wirkende Kraft (Scherkraft) ER 
% Ba pller + Mer {MN rn... nn. (136), 
für die normal zum Rand wirkende Kraft R 
N enpe,.(brr + Mer + MN Le. (137) 
und für die Querkraft 
ERBETEN ne ee 2. SDR 


Für die kartesischen Komponenten des Momentes der im Rand wirkenden Spannungen in 
bezug auf einen beliebigen Raumpunkt mit den Koordinaten a; folgt 


- 


Yu armen SR —aN Sarcn (1 +2wH +wK)dw ... . (139), 


| = 


11) Vgl. H. Neuber, Z. angew. Math. Mech. 23 (1943) S. 323 und 324. 
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bzw. nach Integration unter Beachtung der Gl. (20), (45), (70), (92), (94) 
ein mit A?/12 behafteter Rest vernachlässigt wird: \ 


Hierbei liefert‘das rechte Glied der rechten Seite offenbar das Moment M, in bezug auf 
auf der Schalenmittelfläche befindlichen Randpunkt. Wegen Gl. (131) kann dies in der kürzeren 
Form 


Es folgen ö | ui, 
: B= rk M; = 99 MP rt: a far TE (142), 


Te Mehr 


Bei vorgeschriebenen Werten für 8, N, Q, Bund T würden sich insgesamt 5 Bedingungen er- 
geben. Wie schon Thomsen und Tait!) für die Platte gezeigt haben, ist jedoch die 
Einführung von Ersatzkräften zulässig, wofern diese einem statisch gleichwertigen Kräfte- 
system entsprechen, das durch Zerlegung des Torsionsmomentes in ein Kräftepaar gewonnen 
wird. Zur Ableitung der entsprechenden Beziehungen sei das Moment der Randspannungen 
für ein endliches Randstück näher untersucht. Für das zwischen zwei Randpunkten I und II 
übertragene Moment folgt aus Gl. (140), (142) und (143) 


I 1 
ji ds = fe, rn (9 a PR Be Du Ae re (144). 


Nach Einsetzen von n; aus Gl]. (131) folgt mit dx statt r! ds der Ausdruck 


II II 
ja ds = sm #0 Pr Brida, Du ar aa 0 (145), 


bzw. nach partieller Integration 


II II f 
[ds Ham 2 — ah) [PP N" Mar] + Bri}ds + sum! — a)" |, (146). 


Da zwei Kräftesysteme einander statisch gleichwertig sind, wenn das Moment in «bezug auf 
jede beliebige im Raum befindliche Achse übereinstimmt, so zeigt Gl. (146), daß außer dem 
Biegemoment B im wesentlichen nur die Ersatzkraft j 


APR (NE Pr pe) N (147) 


für die Randbelastung maßgeblich ist. Das noch verbleibende Restglied &.m nn N” Tr 
entspricht zwei senkrecht zur Schale gerichteten Einzelkräften in den Punkten I und II, welche 
sich mit den übrigen äußeren Kräften der Schale zusammenfassen lassen. Ihr Einfluß auf die 
Spannungsverteilung ist nach dem Prinzip von St. Venant vernachlässigbar. Die Kom- 
ponenten der Ersatzkraft ergeben sich mit Bezug auf Gl. (56) zu: | 


IS=S+urrT, | 
* 
N=N Hub Tr» re use (148). 
ei 
Q == Q FR; Te di | 
Zusammen mit dem Biegemoment sind daher nur insgesamt 4 Größen bei Aufstellung der Rand- 
bedingungen zu berücksichtigen. Dies steht im Einklang mit den in der Theorie auftretenden 
Differentialgleichungen, deren Lösungen bei 4 Randbedingungen gerade eindeutig festgelegt 


sind. Entsprechendes gilt, wenn längs des Randes Formänderungs- und Spannungsgrößen 
gleichzeitig vorgeschrieben sind (gemischtes Randwertproblem) 


. 


12)'Handb. d. Tiheoret. Physik, Bd. 1, T.2, Art. 646 ff, 
Eingegangen: 11. Juni 1948. 
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geschrieben werden. Seine Komponenten sind das Biegemoment B und das Torsionsmoment T. 
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| dynamics is geometrically representable by two parts of surfaces each bordered by a closed curve. Four types _ 


. courbe close. On peut distinguer quatre iypes de positions limite, e’esi-ü-dire il y a huit courbes limite qui 


et W. Wwest traitant le m&me sujet sont complöt&s ei corriges. 


gegabelten Verdichtungsstoßes fand kürzlich eine allgemein zugängliche Darstellung durch 
Wuest[4]. Jeder der drei genannten Autoren benutzte neben der „Herzkurvenmethode“ als 
wesentliches Hilfsmittel zur Gewinnung eines Überblickes über die Gesamtheit der möglichen 
Stoßwellengabeln die Aufsuchung ausgearteter Grenzlagen, die sich analytisch durch Grenzüber- 
gänge ermitteln ließen. .Die von Weise und Eggink gewonnenen Ergebnisse wurden durch 
Wuest [3] [4] erweitert und z. T. abgeändert; doch ist auch in [4] der Sachverhalt unvollständig 
und teilweise unrichtig wiedergegeben. i Be: 


nungs- bzw. Konstruktionsverfahren angegeben; es zeigt.sich, daß die 
Grenzlagen sich zur lückenlosen Umrandung der zwei Lösungsbereiche Bild 1. Stoßwellengabel 
zusammenfügen. Die genauere numerische Berechnung machte außerdem 
Einzelheiten erkennbar, die den bisherigen Darstellungen entgingen und a Wirbelschicht 
zu Unstimmigkeiten Anlaß gegeben haben. 


WED ji N 


ir KErVR, vH 


PER % = R | SR a % @ N I N au I z } Et G & 
 . Grenzlagen gegabelter Verdichtungsst 
Von Franz Wecken in Haltingen (Kr. Lörrach) 


Die Gesamtheit der gasdynamisch möglichen stationären Stoßwellengabeln läßt sich geometrisch dar- _ 
‚stellen durch zwei Flächenstücke, die von je einem geschlossenen Linienzuge berandet werden. Es werden 
‚vier Typen von Grenzlagen, also insgesamt 8 Randkurven unterschieden, die an 8 Eckpunkten zusammen- 
stoßen; Berechnungs- bzw. Konstruktionsverfahren und numerische Werte für sämtliche Grenzlagen werden 
angegeben. Frühere Darstellungen von A. Weise, H. Eggink und W. W west werden ergänzt bzw. 
berichtigt. ; * / ah j 


The whole of the stationary bifurcations of shock waves that are possible according to the laws of gas 


of boundary positions are distinguished, that is 8 boundary curves meeting in 8 corners. Methods for compu- 

tation and construction and numerical values for all boundary positions are given. Former treatises by A. 
Weise, H. Eggink, and W. Wwest are compleied and adjusted. / FR 

L’ensemble des bifurcations stationnaires d’ondes de choc &lant possible selon les lois de la dynamique 

. des gaz peut ötre represent& geometriquemeni par deux parties de surface dont chacune est limitee par une 


se rencontrent en huit coins. L’auteur publie des methodes pour le calcul et la construction. De plus il donne 
des valeurs num£riques pour toutes les positions limite. Lestravauz anterieurespar A.Weise, H.Egg i nk 


Mno:kecTBo Bcex CTAIHOHAPHEIX PA3BeTBIeHnÄa BONH OT YAapa, HONYUHAIOIIMXCH BaAKOHaAM 
'Ta580B0f AHHAMHKH, MOKeT ÖBITB TEOMETPHYECKH IIPeNCTaBlIeHO ABYMA HOBEPXHOCTAMH, 
KaKad US KOTOPEIX OFPaHuyeHa HEKOTOpoü 3amkHyroi kpusoli. Pasımuamrea yersipe 
TaNa TPaAHHYHEIX IIOJIORKeHMÜ, T. €. BCerO 8 TPaHHYHLX KPHBEIX, CXOMMMHXCH B BOCBEMM 
yrıax. ]lamrca pacyeTHke MIN TeOMeTPHUYecKHe METONBI MH YHCHIEHHBIe BEIIMYHHE IA BCEX - 
TPaHHYHBIX MONOoKeHnü. JlonosHamrca HAIM HCHPaBıamTcaA NnpeikHmne pabors A. Beüse, 
X. Iıruuk u B. Boer. 2 RE re 


| $1. Einleitung . ” 
-Die von Weise [l1]und Eggink [2] unabhängig voneinander begründete Theorie des 


Wir unterscheiden im folgenden vier mögliche Typen von Grenzlagen (vgl. $2); von diesen . "ER 


sind die ersten zwei (vgl. $3) schon in [1] und [2] grundsätzlich richtig angegeben. Der dritte 
und vierte Typ wurden nirgends klar getrennt (freilich kann man ihre‘ 
Verschiedenheit aus einigen Figuren herauslesen); Formeln wurden nur 
für den dritten Typ (&=0, $4) in [1] aufgestellt und erlaubten nur 
Auflösung durch Probieren. — In [3] und [4] wurde ein fünfter Typ 
fälschlich als vorhanden behauptet. 


In dieser Arbeit werden für alle vier Typen handliche Berech- 


h Hauptstoß, a und 
b Nebenstöße, 


In den Bezeichnungen folgen wir z. T.Wuest[4] (vgl. Bild 1). Die Gebiete 3 und 4 stimmen 


hinsichtlich Druck und Strömungsrichtung überein. Die drei Stoßfronten a, b, h kennzeichnen 
wir durch ihre Druckverhältnisse 


bzw. deren Logarithmen 


A=—ın&, A=—In&, Ah=— In &: 


es ist 
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| schwindigkeitsmaß benutzen wir neben der Manu Zahl a; dien mit der k 
Ba gebildete Zahl M*; es ist 


Die zwischen M,, M, und &, bestehende Beziehung. 
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mM; BvR | 
läßt sich nach jeder 3 drei Größen ee auflösen; wir können z.B. schreiben 
M;=Y(M,A) :- 2.200 nn Sa 


Entsprechendes gilt für die Stöße b und h. — Stoßwinkel y und Umlenkung ö (Bild 1) schreiben 
sich in Funktion von ‘Anströmung M und Druckverhältnis & (bzw. von M und a 


A A) 3 a 
y=aresin ( Ne = EL E <. : 
k ü 
+ Var ar! H 
s=artg —— —— = Q(M.A).. .....- (4). 
u } Mr 6 


(4) ist bei festem M die Gleichung der ‚‚Herzkurve“, 


$2. Die Fläche G | 
Ein Wertesystem M,, &., & bestimmt dann und nur dann eine Stoßgabel, wenn es die 
Gleichung ö, + 6, = ö„oder 
Q(M,„A) +2(PM,A),A)= RM, A+A)- -: 0... . (5) 


und die Ungleichungen.(1) erfüllt, (Die physikalisch ebenfalls wesentlichen Ungleichungen M, > 1, | 
M, > 1 folgen mittels (2) und (4) aus (1), brauchen also nicht aufgeführt zu werden.) (5) wird im 3 


1 4 
&u- & Raum durch eine reelle algebraische!) Fläche dargestellt; aus dieser wird durch i 
1 


(1) ein ganz im Endlichen liegender Teil @ herausgeschnitten, der evtl. aus mehreren Stücken be- 

stehen kann. Die sämtlichen möglichen Gabeln entsprechen eineindeutig den Punkten von @. ö 
Als Berandung von @ haben wir nach (1) vier Arten von Kurven zu erwarten, je nachdem 

&, oder &, gleich Null oder Eins wird, sowie vier Arten von Eckpunkten: 


Aal B: &=1, &=0 | : 
0: == DYn mn, Be 1 F 9. 


Die Herzkurvenmethode nach Weise und Eggink (vgl. Bild 2in [4]) liefert graphisch 
alle Lösungen A, von (5) zu gegebenen M, und Aa; die Zahl n = n (M,, A.) dieser Lösungen gibt . 
zugleich die Vielfachheit an, mit der ein Punkt der M,-£a- oder M,-M,-Ebene durch die Pro- 
jektion von @ auf diese Ebene überlagert wird. Es liegt nah, @ in einem Alle beiden Diagramme 
darzustellen; beide Abbildungen sind, ohne Ränder, topologisch äquivalent wegen der Umkehrbar- 
keit von (2). Wenn jedoch (was sich unten zeigen wird) alle vier Randkurven existieren, so sind 
diese Diagramme zur Veranschaulichung von @ unzweckmäßig, da offenbar stets eine Rand- 


kurve auf einen Punkt zusammenschrumpft: im ersten Falle der Rand Ben 0 SR den Punkt 
1 ‚ Bi: 

ee n, = ns im zweiten Falle (Bild 13) &,=0 auf den Punkt = M ( oder 

M* — M*— 7): 


Yk 


‘) Obwohl (5) in transzendenter Form geschrieben ist! 
J 


An dieser Stelle wird die Zuordnung der Gabeln zu den Dingrann 


EEE N Melde Man man Sul 


PH: een Parametern M,, M, aus. Aus der graphischen ‚Konstruktion der Lösungen von (5) a 
BE. ‚405 erschließt man ahnlich folgende Tatsachen: 


Mm, und M, der ne (d. h. ee mit kleinstem n vom zweiten stets um ein endliches Stuck 


? menhängenden Gebiet der M,-M,-Ebene. ER 


"destens) zwei getrennten Teilen G, und @, besteht; G, ist einfach zusam- | 
_ menhängend und projiziert sich ehne Raid eineindeutig auf ein Gebiet 


. zuG, gehörenden Lösungen, so daßB@=G, +G, ist. Über die Gestalt von Herzkurven 
-@, sei zunächst nichts ausgesagt. 
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'hauptet wird. Einen klareren Überblick über die Struktur von G, und seine Darstellung in den 


1er sp äter ($ a ein Bi 


den een sr zur Klarinz des er, von G Aunächen) von 1 den za 


EA Liegen 2 Schnittpunkte der Herzkurven vor (Bild 2), so bleibt bei stetiger N Br 


entfernt. . 13227 } 
-b). Die erste Var ist nalen und stetig in einem einfach zusam- 


Hieraus folgt (was schon in [4] hervorgehoben ae, daß @ aus (min-. 


der M,-M,-Ebene. Mit @, bezeichnen wir vorläufig die Gesamtheit der nicht Re 


ES Br: 


Da der zweite Schnittpunkt nur zusammen mit dem ersten ee kann, erscheint im f4 ce 
M,-M,-Diagramm G, als Teilbereich von @, (vgl. Bild 13). Daß @, und @, getrennt liegen, ird 
erst anschaulich eient z. B. bei Benutzung der Koordinaten MZ und 6. Egginkf[2] zeichnete 
Kurvenscharen M, = const in ein M}-6,-Polardiagramm, wobei zwei getrennt liegende Scharen 
(entsprechend G, und @,) erkennbar sind, aber wegen der z. T. fehlenden Ränder die genaue Er- 
streckung von 6, und @, nicht klar wird. Wuest übernahm das gleiche Diagramm ([4], Bild 6) 
und zeichnete A Ränder ein, dagegen keine Kurvenscharen; erst beide Darstellungen zusammen- 
genommen machen das Gesagte deutlich. 


Hinsichtlich @, glaubte man das gleiche wie für G, zu finden: schlichte Abbildung auf ein 
einfach zusammenhängendes Gebiet der M,-M,-Ebene (folglich einfacher Zusammenhang bei 
einfach geschlossenem Rand). Doch ist in anchen Fällen selbst bei guter Zeichengenauigkeit die 
genaue Lage und selbst die Zahl n der Herzkurvenschnittpunkte schwer zu ermitteln, und wir 
werden in der Tat finden, daß es auch dreischnittige Lagen gibt, während in [1], [2], [4]n <2be- 


verschiedenen Diagrammen, als ihn die bisherigen Methoden ergaben, erhält man aus einer voll-. 5 
ständigen Untersuchung des Randes. Kennt man dessen genaue EN so ar es leicht, in ihn as 
Flächenstück einzuspannen. ' 


Zur graphischen Veranschaulichung = man möglichst ein Diagramm ablen) in dem das 
fragliche Flächenstück mit Rand topologisch abgebildet erscheint. Das gilt für @, von keinem 
der bisher benutzten Diagramme, da in diesen nicht einmal der Rand von @, ee ab- 
gebildet ist. Nimmt man £&, und £&, als Koordinaten, so ist zunächst klar, daß @, sowie @, ins i 
Innere des durch (1) gegebenen Quadrates Q abgebildet werden, und die Ränder von @, und G, 


‘auf den doppelpunktfreien Rand von Q. Man kann daher erwarten, daß die Abbildungen 5 >Q  , 


und @,—0@ (Bild 14, 15) mit Rand topologisch sind, was sich in der Tat bestätigt (vgl. $6). Zum 
Nachweis haben wir zunächst die Randabbildung zu untersuchen, was wir für G@, und @, gemeinsam 
durchführen. 


$ 3. Die Grenzfälle &=1 und 4 =1 
Bei &=1,4.=0 („vordere Grenzlage“ nach Wuest) wird 2(M,„4A.)=0 
Y(M,,A.) = M,; aus (5) wird daher durch Grenzübergang A.—0: 


02 02 oy 


| 20 
Ur) + MA MOM A) (7). 


‚Dies ist eine implizite Gleichung zwischen M, = M, und A, =4, oder & = £&, die in anderer 


Schreibweise auch in [1], [2], [3], [4] angegeben ist?). Ihre explizite Auflösung dürfte unmöglich 
sein; die numerische Ermittelung der Lösungskurven ist mühsam, wird aber dadurch erleichtert, 
daß man ihren ungefähren Verlauf aus dem Herzkurvenverfahren erschließen kann und außerdem 


2) Die Formeln (13) und (14) in [14] enthalten mehrere Fehler. 


% I 7 [ 
er ni 
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Ei t Rat „ee . 
F en and ve En f 
r h ur "RR Ki: ie) un 8 n“ B 


RZ HR, 


|  Wecken, { 


« 6- x , 2 


e \ Bild 3 darstellt. Aus zeichentechnischen Gründen sind &, und M? as 


B, und B, fallen dabei zusammen. nr Bad, 


i - A | 
Bi N: ider | n wächst &, monoton von O bis 1,d.h A, B, und A, B, je für nd ne 
Br se a =, abgebildet. My hingegen ändert sichnurauf A, Bamonoton 
En und erreicht auf A, B, ein Minimum in E. 2 Ba, 
r N "Als eindeutige Funktionen von &, auf A, B, und A, B; erhält man nach (3), (4) die weiteren x R 4 
Er Gabelparameter y„ = y, und 6, = ö,; weiter ist y„ = are sın (1/M)). Es ist nur zu beachten, S: 
Br i daß ö, und 5 — y„ auf SB, negativ, sonst p- 


„ sitiv zu wählen sind. In S berührt A, B; die 


as Kurve 2(M,Aı) = 0 oder 
1 
“ PFERBR ed. 
Y; | ı 1+6°’. 
| u” IA FRE: 2 j 
Bid 3. Randkurven &, =1. Bild 4. die in’ Bild 3 z. T. eingezeichnet ist. 
hy PEN DE ‘ Numerische Werte für die benannten 


RR ST Punkte sowie die entsprechenden Bezeichnun- 


gen in [4]finden sich in Tafel 2, Formeln in Tafel 1. 
Bei A,—0 (‚hintere Grenzlage‘‘) erhält man aus (5) 


202 02 
En (7 (M,,4.),0) = a Mr Au) era seer (8) 
1% als implizite Gleichung zwischen M, und A„—= Az. Es ist in diesem Falle auch eine explizite Dar- 
stellung möglich: ee 

Eu 1 a) 
I er f lebe | | 

M2— a 1 a . — e +1 (9) 
ae +1, M=75 EIER 
mit d| 
2k 8a Hi 
FE TE ET ET 2 N 
1+5+ Ya tee (147) { 
das obere bzw. untere Vorzeichen gilt für @, bzw. @,. Daß (8) durch (9) identisch befriedigt wird, | 


verifiziert man durch Einsetzen. Nach (9) sind M, und M, und damit auch die übrigen Gabel- N 
parameter eindeutige stetige Funktionen von & (0 < &. < 1) für jeden der beiden Lösungsbereiche; | 
damit haben wir die gewünschte Darstellung der Randkurven A, D, und A,D,. In Bild 4 ist 
für beide Ränder M7? über &, aufgetragen; dabei fallen hier D, und D, zusammen. Wieder ver- \ 
läuft M, monoton auf A, D,, erreicht aber ein Minimum (F) auf 4,D,. M,ist auf beiden Rand- 

kurven monoton. A, und A, sind die uns schon bekannten Eckpunkte, 


$4. Die Grenztälle &&=0 und ,—0 
Während für die Ränder &,=1 und &,=1 die Herzkurvenmethode versagte und ana- 
Iytische Verfahren Anwendung fanden, ist sie für die bisher weniger beachteten Ränder =0 


und &,= 0 mit einer geringen Modifikation wieder verwendbar. Dazu brauchen wir die Grenz- 
lage der Herzkurve für M), = «. 


Während die Stoßpolare für M— c bekanntlich in einen Kreis übergeht, bekommt man 
für die Herzkurve bei üblicher Darstellung keine Grenzlage, da bei M),—> ® für endliches A immer 
ö= 0 wird. Die Herzkurven schmiegen sich bei MÄ— ® mehr und mehr der positiven A-Achse 


an (Bild 5). Verschiebt man jedoch die Kurven vertikal so, daß die Scheitelpunkte zusammen- 
fallen (Bild 6), indem man über ö anstatt A eine Größe 


ae! ae?) ERRER ER =In (M: (1+ k) Rn k) 


ä #4 en Bild 5. Se ns ? : 4 Bild 6. Herzkurvenschar mit R - Bild. ?. Herzkurvenschnitte a Ss 
er; mit gemeinsamem Fußpunkt gemeinsamem Scheitelpunkt - } bei &, =0 : SE 
E yr Grenzherzkurve H, für Mı =» Zweischnittige Lage Se 


Beide Tatsachen lassen sich en: INC Dreh Die Gleichung von H, geben wirin 


zwei "Gestalten: 


Be a Ba Re BL 3 
ee 4% in x sin bi ee Are ge EN EN ' 


Stoßwinkel und u Raunng, sind gegeben durch. 


sin?y= en Mm = a Rn RR 


k TER Se 
Legt man das Herzkurvendiagramm Nach Bild 6 zugrunde, wobei H,„ nach (10) berechnet ist, B 
so läßt sich das Herzkurvenverfahren sofort auf den Fall M. [> &, M,< © ausdehnen (Bild 7). ° 


# 


x 


„0a, Ös, 6) und A, sind aus der Figur abzulesen; y. und yz ergeben sich nach (10'”), 9 nach 8). a 
Es gibt keinen, einen oder zwei Schnittpunkte, je nach der Lage von N bzw. dem Werte des Para- 


Funktionen von A’ oder von M, erscheinen. 
Trägt man beide Kurven in ein &-y„-System ein (Bild 8), so findet man, daß auf jeder von 
"ihnen &, monoton von 0 bis 1 läuft. Die Endpunkte D,, D, sind die schon bekannten. Die den 
Endpunkten entsprechenden Lagen von N (Bild 7) sind: für C, und 0, M’—=— 9; für.D, der? 
Punkt der größten Ablenkung auf H,; für D,, dazwischen liegend: renze der Ein- und Zwei- 


EN a N hl ne win ibn. Aa 1 


TEEN 


> 


und ö, = 0 ergibt.) 


$ = 05 7 


Bild 9. Herzkurvenschnitte Bild 10. Bild 11. 


Bild 8. Randkurven E,=0. : \ 
a T senkrechter Haupt- DER See LEN ED" Grenzfall & = 0, Grenzfall &, = 0, 
B stoß great BuyparZyre, Eymız erste Lösung. &, = 0,5 zweite Lösung. &, = 0,5 
metrische Schnittpunkte 


Um die noch fehlenden Ränder B; 0; oder &= 0 aufzusuchen, haben wir offenbar für 
Haupt- und Nebenherzkurve die Grenzform H, zu nehmen, denn es ist M=-&M = % 
Für endliches &, ist A, = —o, daher 6, = 0 und „= 0. Die beiden kongruenten Kurven H,„ 
und Ha sind also vertikal gegeneinander zu verschieben (Bild 9). Es gibt stets zwei Schnitt- 
punkte in symmetrischer Lage, an denen man ö (= 6,), A, und Ay abliest; nach (10”) folgen 
y» und yr. Die Verschiebüng d, beträgt nicht, wie man vermuten könnte, A.. Bezeichnet man ' 


2 3) Die Gleichungen (10) kommen in anderer Gestalt schon in [1] vor. 


Pr = xsin?ö AT + cos öy I — es : sen (10°). | 2 


\ \ J rg = 
{ 1 € —4’ 1 : Te 
4 # ar -) 


Hier ist offenbar &, = 0, &, #0. Der Fußpunkt N der Nebenherzkurve wird durch (10°) gegeben. 8 | 


meters A’ für N. Man erhält also zwei Randkurven, auf denen alle Gabelgrößen als eindeutige E7 


schnittigkeit. (Formeln und Zahlenwerte vgl. $5, auch für den Punkt T, der sich aus & = 0 ER. 5 


Kr \ Jl "R a RSEH ni 


Ei 0 ri. 2, &,) die vertikale Verschiebung der Schieiunkte von PR? R a i urve 

mr bei endlichem M,, so wird nach einfacher POSRIRE u a | Hi 

Be TE: ER 2 

R >; 

e: 2 R d, DEAL N +8 - leid). Asaalıe a 

: # ‘Man en also &, als unabhängigen Para wählen, hieraus nach aı d, und damit Bi “s u | 
ES, ' Bild 9 die sämtlichen Größen auf B, C,und B, C,in Funktion von &,ermitteln. Auch die Intervall- | 
Be enden B; und C; ergeben sich unmittelbar, ohne neuen Grenzübergang. Die beiden für irgendein an 
0% &ugefundenen Lösungen (Bild 10 und 11) unterscheiden sich nur durch die Vorzeichen von = &, 
LAN #—yn und $—y,. In beiden Fällen herrscht in den Gebieten 1 und 2 der Druck Null, dehhat 
D 'Pılpz den bestimmten Wert &,. Die unendlich starken Stöße k und b ergeben gleiche Umlenkung 

Er bei verschiedenem Stoßwinkel. — Das graphische Verfahren (Bild 9) ist übrigens nützlich für die 


ui - 5 Anschauung, aber entbehrlich; nimmt man z.B. ö, als Parameter, so erhält man aus (10°) die 
2.0 beiden Werte A; und A;, aus (10°) y, und y, und aus (11) &,, denn es ist offenbar d, = Ay, — 4}. 
N Wegen dieser einfachen Berechnungsmöglichkeit verzichten wir für die Ränder &, = 0 auf graphi- 
sche Darstellungen. Alle erwähnten Parameter laufen auf B,C, und auf B,C, monoton. 


$ 5. Die Eckpunkte 
Daß die aufgefundenen zweimal vier Randkurven den gesamten Rand von @ bilden, folgt 
daraus, daß sie sich zu zwei geschlossenen stetigen Kurvenzügen 4A, B,0,D,A, und 4, B, 6, D,As 
zusammenfügen. Die stetige Abhängigkeit der Gabelparameter von &, bzw. &, auf den einzelnen 
Bögen A; B; usw. haben wir bereits eingesehen; es bleibt noch zu zeigen, daß sie auch an den 8 


Tafel1. Eckpunkte von G und Randpunkte mit ö, = 0 
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sin du = ein = ap, = Vr 2“ 
1 

& mE sen 

sh ART 


Sin da = Ge + Dane 


ANA 1 1 1 1,245 | 1,2 0 0 0 53,46 | 53,46 

AR U 1 1 1 2,565 | 2,565 0 0 0 22,95 | 22,95 

b, E 1 0 0 oo oo 0 —45,38 | —45,38 0 112,3 

DZ 1 0 0 [e0) 00 0 45,38 45,38 0 67,69 

Q, E' 0) 0) 0 oo to 0 —19,84 | —19,84 0) 155,8 94,32 
(, W 0 0 N) 00 co 0) 19,84 19,84 0 24,15 | 85,68 
Di’+@ 0 1 0 (6) 1 45,38 0 45,38 | 67,69 | 90 67,69 
Ds Y 0 1 0 00 1,471 | 42,81 0 42,81 | 57,76.| 42,81 | 57,76 
S} D l 0,416 | 0,416 | 1,484 | 1,484 0 0 0) 42,36 | 42,36 | 90 

Ah F 0 0,139 0 00 5,632 | 18,00 | *—18,00 0 21,86 | 153,5 90 

E — 1 0,5 0,5 2,097 | 2,097 0 19,4 19,4 28,5 40,5 40,5 
F — | 0,65 1 0,65 2,480 | 2,198 | 6,9 0 6,9 29,2 26,9 29,2 
H C 1 0,583 | 0,583 | 1,334 |. 1,334 0 7,39 7 


‘) Spalte 1: Bezeichnungen dieser Arbeit; Spalte 2: Bezeichnungen nach Wuest B]) 4]. 


| | ı an 
Mr zur en, de Gabeln üblichen Parameter behaupten, wir die Steige 
Die Tatsache, daß m. a. W. von ‚den 16 End- 
punkten der 8 Kurven je 2 identisch sind, haben 
WEL, durch die Bezeichnungsweise der Punkte 
schon vorweggenommen. "Statt eines besonderen Br 
 Nachweises stellen wir in obiger Tafeli die 
A eraur, Berechnung der Eckpunkte sowie der Punkte 
8483) und 7 ($4) erforderlichen Formeln zusam- 
er men, die großenteils aus [1] und [4] übernommen 
Br sind. Das obere Vorzeichen gilt jeweils für @,, das 
untere für G,. Bi LundSity=yn=ypn=y 
2.2 bzw. M\= M,—= M gesetzt. — Tafel 2 zeigt Zah- | - } 
lenwerte für «= 1,405, die nach Tafel 1 berechnet pua 12. Zweitach entartete Grenztälle (Hekpunkte). Die 
e' sind. Hier sind auch Werte für weitere Punkte EZ, _ Stoßtronten sind z. T. unendlich stark (stark ausgezo- 
-F, H des Randes (vgl. Bild 3, 4, 17) aufgeführt, EEE ER 
‚für die explizite Formeln nicht vorliegen. — Für die 8 zeigt Bild 12 die geometrische ea 
Gestalt der Gabeln. Die Grenzlagen mit &, = 0, sowohl die einfach- wie die zweifach-entarteten, 
| sind natürlich physikalisch nicht realisierbar, man kann ihnen aber beliebig nahe kommen, sofern 
überhaupt die dem Inneren von @ er nicht entarteten Gabeln als reell gelten können. 


. | S 6. ade Nungen. von Gı Ha GN 

B ‚Aus dem vorstehenden folgt, daß die am Ende von $ 2 erwähnte Abbildung von @ auf R: hin- 
sichtlich der Ränder von @, und von @, (je für sich) topologisch ist. Um das gleiche auch für das 
Innere einzusehen, Peehlen wir ebene naher Darstellungen von @i im M,-M;-Diagramm 


. Bild 13) sowie von @, (Bild 14) und @, (Bild 15) im &,-&-Diagramm. 
8 
E In Bild 13 erscheinen, wie schon erwähnt, die Bö- 0 Le: 
A ; gen B;C; als ein Punkt B= (0; 4, D, werden als Kurven, 
Fe A; B; und 0, D; als Teile von Geraden abgebildet. Das Bild 
7 des Randes von @,ist eine Jordan-Kurve A,D, BA; 
i 1 AL s 
0 j 
A a 
/ 2 i 
A KM? D, > 
£ 0m 08 7 | 
; ° Bild 13. %M}?- k M32-Diagramm. «  Bild14. &,- &;-Darstellung von @.. 
2 n = Schnittzahl des betreffenden Mı = const;, ———-— M: = const; 
| ; BPFEEHETERPARAER 1-2: 5, = 0 (ST) 


für G, jedoch ergibt sich ein Linienzug 4,EBD,FA,, in den sich offenbar kein schlichtes 
Flächenstück einspannen läßt. In der Nachbarschaft der Strecke A,E (Bild 13) muß die Be- 
deckungszahl für @, allein mindestens 2, für @=@,+G, daher n2 3 sein (vgl. $2). Die ein- 
H fachste Annahme ist die, daß eine Faltungslinie EX ein Gebiet mit n = 3 begrenzt, wobei die. 
Lage von X auf A,D, noch zu bestimmen ist, und daß @, schlicht, abgebildet ist. 

Es muß also Wertepaare M,, M; geben, für die die Herzkurven drei wesentliche Schnittpunkte 
liefern, und dies wurde in der Tat bestätigt gefunden bei M, = 2,46, M, = 2,22 und 2,01 (vgl. 
Bild 16). Bei Erreichung ‘der Linie EX (Bild 13) vom dreischnittigen Gebiet aus rücken zwei 
Schnittpunkte zu einem Berührungspunkt zusammen. 

Zur genauen Ermittelung von X überlegt man sich leicht: dem Punkt X entspricht 
eine Berührung beider Herzkurven von zweiter Ordnung im Nebenherzfußpunkt; es muß 
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- 154 | Wecken, 
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dort also n ®2 pm, As)» = 77 (Mu Ae) oe N ee a N) u. 
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sein, während die linke Seite von (12) auf 4,X kleiner, auf D,X größer als die rechte ist. Da beide De | 


Seiten von (12) nach (9) auf A,D, as allein von £, sind, läßt sich X aus (12) besikinEne Ang 
Dies wurde numerisch nicht durchgeführt.) | j ni RE 
; Nach Bild 13 liegt die niedrigste Anströmgeschwindigkeit für das zweischnittige Gebiet on 
nicht, wie bisher angegeben, bei M, = 2,565 (d.h. A,), sondern bei M,= 2,480 (d.h. A, und 
für den zweiten Lösungsbereich überhaupt, zugleich für 
das dreischnittige Gebiet, bei M, = 2,097 (d.h. bei EB). 
Um die gewonnene Vorstellung zu erhärten und, 
eine auch für @, schlichte Darstellung zu erhalten 
haben wir die durch Herzkurvenschnitte. gewonnenen 
Gabeln‘), für @, und G, getrennt, in ein &-&-Diagramm 
eingetragen (Bild 14, 15) und durch Kurven M, = const 
und M,= const verbunden®) unter Berücksichtigung 
der in $$3—5 berechneten Randwerte. Obwohl sich aus 


Bild 16. Dreischnittiger 


Bild 15. &,t&, -Darstellung von G.. Kurven Fall (schematisch). Ein 
verschiedener Scharen schneiden sich z.T. Schnittpunkt zu @,, zwei 
doppelt zu G: gehörig 


endlich vielen Punkten und Kurven mathematisch streng die Gestalt’einer Fläche nicht ermitteln 
läßt, wird man aus diesen Diagrammen folgendes erschließen: 

Durch jeden Punkt in Bild 14 ‘und 15 geht eine und nur eine Kurve M, = const sowie eine 
und nur eine M, = const; zu jedem Punkt des Quadrates gibt es daher je eine und nur eine Gabel 
von @, und von @,, oder: G, und @, sind auf Q mit Rand topologisch abgebildet. M.a. W.: auf 
G, und auf @, sind alle Gabelparameter eindeutige stetige Funktionen von &, und &. — Daß 
zwei Kurven M, = const, M, = const sich im Innern von Bild 14 höchstens einmal schneiden, 
ist gleichbedeutend mit der in $2 erwähnten Schlichtheit der Darstellung von @, (ohne Rand) 
in Bild 13. — Ferner: Bild 15 wird durch die Linie EX in zwei Teile geteilt, deren jeder für sich 
‚in Bild 13 schlicht dargestellt ist; dabei entspricht jedem Punkt in dem Dreieck EA, X ein anderer 
Punkt außerhalb davon mit gleichen Werten für M,, M, und &.. Die durch die Punkte von EX 
(in Bild 15) laufenden Kurven M, = const und M,= const haben hier eine gemeinsame und 
zwar vertikale Tangente. — Die Kurven M, = const und M, = const auf @, bestehen z. T. aus 
zwei getrennten Teilen, während sie auf @, stets zusammenhängend sind. — Die von Wuest [4] 
fälschlich als Randkurve des zweiten Bereiches, bezeichnete Linie M, = 1,8417 (kM}: = 0,5712) 
verläuft in Wahrheit (vgl. Bild 15) von & = 1, & = 0,2 durch das Innere von @, nach & = 0,4, 
&=1; dazu kommt noch A, als isolierter Punkt. 

Um irgendeine andere zweidimensionale Darstellung von @ zu erhalten, geht man zweck- 
mäßig von der Vorstellung der beiden Quadrate (Bild 14, 15) aus und bildet zunächst nach $$ 3—5 
die Ränder ab. Um etwaige Faltungslinien und damit die genauen Konturen und Bedeckungs- 
zahlen zu bestimmen, muß man eine Kürvenschar, z.B. M, = const, zu Hilfe nehmen. Bild 17 
und 18 zeigen die so sich ergebenden Darstellungen in den Koordinaten 65, 6. und M*, Yy„ unter 
Fortlassung der Kurvenscharen; die entsprechenden Figuren in [4] (Bild 3 und 5) sind hinsichtlich 
G@, fehlerhaft. Bild 17 liegt im wesentlichen schon bei Weise ([1], Bild 10) vor. Punkt H, also 
das Maximum von ö, auf A,8, ist in Tafel 2 aufgeführt. V wurde nicht berechnet. 

In [5] wurde die Theorie des gegabelten Verdichtungsstoßes auf den Mach- Effekt an- 
gewendet und dabei die Ergebnisse dieser Arbeit z. T. ohne Begründung skizziert. Dort fand 


(Bild 7) ein Diagramm Verwendung, in dem neben &, der Winkel & = ya +3 — yy als Koordinate 
diente, 


*) Entsprechend einem quadratischen Punktgitter in Bild 13 mit der Maschenweite 0,05. Die Zahlen 
an den Kurven in Bild 14 und 15 sind die Werte von 20%k Mt! bzw. 20 k.M%*. 
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Bild 17. Ränder, Faltungslinien. und 

Bedeckungszahlen im ö7, - ö,„-Diagramm 
(Nach Wuest, berichtigt) 
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5 Bild18. Konturen im M} -yy-Diagramm 
g . . Faltungslinie 7 C, geschätzt 


(Nach Wues t, berichtigt) 


zeigen z.B. die Formeln für A, in Tafel 1. Diese werden für «> & sinnlos; für <= = rd 


B und € zusammenfällt, in Bild 18 in die rechte untere Ecke des bis y, 0 fortgesetzten 2 
. Rechtecks rückt. | RR 


3. 


k= T’ für den Punkt 4, also „= y=ym=0, M,=M,=, so daß A, in Bild 13 mit 
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KLEINE MITTEILUNGEN 


Über Integrale von Funktionen, die Produkte 
BESSELscher Funktionen enthalten. 


Jeder beliebige, sich als Welle ausbreitende Schwin- y 


gungsvorgang der Optik, der Akustik, der _Elektro- 
magnetik oder der Wellenmechanik läßt sich bekannt- 
lich durch einen Integralausdruck der Form 


ivt ik le. 
ae! 


worin 
I.» .] Sil@pceosd+yp sin #cos +2psindsin 9—f(9, 9)] 


sinadddp 


darstellen !), falls es sich um „monochromatische“ 
Schwingungen der Frequenz » handelt, oder allgemei- 
ner: durch einen Integralausdruck, der sich aus dem 
Vorstehenden durch Multiplikation mit einer weiteren 
Funktion x (v) und Integration über » ergibt. Hierbei 
sind f(®,p) und y (9, p) geeignet zu wählende Funk- 
tionen, die — z.B. in der Wellenoptik — ganz be- 
stimmt-angebbare physikalische Bedeutung besitzen. 


Bei der Auswertung dieser und ähnlicher Integral- 
ausdrücke der theoretischen Physik ist es oft notwendig, 
eine Reihenentwicklung anzuwenden. Man kommt so 
in vielen Fällen auf Reihen Besselscher Funktionen 
wachsender Ordnungszahl, in denen diese Funktionen 
linear auftreten. Handelt es sich nun um die Bestim- 
mung der Energie oder der Intensität des betreffenden 
Schwingungsvorganges, so erhält man durch Multipli- 


1) Siehe z.B. J. Picht, Ann. d. Phys. (4) 77, 785 (685) 
1925. [Diese Arbeit enthält auch — 8. 752—755 — Tabellen 


82" 79 1(8) s2rt1 7911410) 
n* 11 (R+1)N : 
(n* 11 = 2.4.6.. .2n), (nt = 1.3.5... (2n—1)), auf die hier hin- 


gewiesen sei, da sie im ‚Verzeichnis berechneter Funktionen- 
tafeln‘“ nicht angegeben sind]. 


der Funktionswerte 


kation mit dem konjugiert komplexen Werte jener 
Reihe wieder eine Reihe, deren Glieder in den Bessel- 
schen Funktionen J, bilinear sind. 

Beschränkt man sich nieht auf die Bestimmung der 
Energie (Intensität) in den einzelnen Punkten des Rau- 
mes, sondern fragt etwa nach der Gesamtenergie, die 


(in jedem Augenblick oder im zeitlichen Mittel) durch 


eine gegebene Fläche, etwa durch einen Teil einer Ebene 
hindurchströmt, so ist: eine Integration über das Pro- 
dukt der vorstehend angedeuteteri Reihen von in J„ 
bilinearen Ausdrücken mit dem Flächenelement do 
(= ode dp) auszuführen. Es treten so in vielen Fällen 
Integrale der Form r ; 


Kg 
x. 


auf. Integrale dieser Art: lassen sich nun häufig durch 
wiederholte partielle Integrationen ‘auswerten, was 
sicher auch schon oft durchgeführt wurde. Mir ist aber 
bisher keine Literaturstelle bekannt geworden, an denen 
eine Zusammenstellung der Integrale oben angegebener 
Form zu finden ist. Auch in ‚„Jahnke-Emde“ sind diese 
Integrale nicht aufgeführt. Es sei deshalb gestattet, 
hier die zugehörigen allgemeinen Formeln mitzuteilen, 
ohne daß die z. T. etwas umständliche Herleitung des 
Ergebnisses in allen Fällen mit angegeben werde. Die 
Gültigkeit der mitgeteilten Beziehungen ergibt sich 
unter Berücksichtigung der bekannten Formeln für die 
Differentiation der Besselschen Funktionen 


ddp 
da. 
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wobei n# 0 und n #1 sein muß. 
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Die vorstehende Formel geht für EN m in die vor- Wir geben nachstehend noch — u. zw. mit Ableit 
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Durch VER Anwendung vorstehender Bezie- 


Jm JInJm Jnd 
‚non - Br rs ... ergeben sich fol- 


. gende Beziehungen, wenn wir noch im Interesse der 
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Durch Addition erhalten EN so 
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Die in (2) rechts ei Summe bricht nun 
leider nicht von selbst ab. Es ist daher notwendig, die 
Summe nur bis » = qg—1 zu führen und das Restglied 
in anderer Weise umzuformen. Wir erhalten so 


; InJm d IS: 1 er mt 
Br Ste - 7v 
= da 77 (n+m+g+1-2u) ° 
-1 
1. De (q-1) 
2 Se ee 
i ZI (ntm+g+1- 24) 
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worin wieder 


i . ; IE a 13 Inzgqt1Jm-2. 
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AN Sr (=; BR In-q+1+4Jm-1 
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ng (m-A)2-(n-g+1+A)? (Jn-gq+1+2 Jm-1-1 
— In-g+ Tu) i 


Hier können wir noch den Nenner des zweiten Gliedes 
der rechten Seite umformen, und zwar ist 
m—(n—g+ 1° —2(m+n—g+1)A 

= (m+n—g+1) (m—-n+g—1— 22), 


sich ergibt: 
ar ! : 
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Kritische Betrachtiing des dynamischen Wi- 
„geteilten Platte bei Stoßbeanspruchung. 
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Potsdam. Babelsberg 2. Johannes Aue 


derstandes einer in mehrere Schichten auf- 


I. Die dynamische Beanspruchung 
einer Platte. 


Eine Platte werde stoßärtig beansprucht: 


a) durch große Massen geringer Geschwindigkeit, 
b) durch kleine Massen hoher a 


In jedem Fallist die Beanspruchungeine d. amische, 
die kinetische Energie oder Wucht einer auftreffenden 
Masse m der Geschwindigkeit v, ist 


; 
gege ..() 


Wenn die Masse senkrecht auf die Platte auttrifftund - 
kurze Zeit später einen Moment zur Ruhe gekommen 2 “E 
ist, hat in diesem Moment die Platte die gesamte Ener- PER 
gie bzw. Arbeitsfähigkeit der Stoßmasse aufgenommen. 

Gefordert wird, daß die Platte diese Energie ohne 
Bruch aufnehmen kann. 

Die betrachtete Stoßmasse m drückt während ihres 
Auftreffens mit einer Kraft P auf die Platte, die pro- 
portional der Geschwindigkeitsverzögerung der Stoß- 
masse ist 


und die an die Platte abgegebene Energie bzw. Arbeit 

ist die Summe über die ir jedem Moment der Berüh- 

rung wirkenden Stoßkraft maldem Durchbiegungsweg s PR 
der Platte 


A=/[P(s)-ds 


P(s) ist die Druckkraft als Funktion vom Biegeweg s, 
der seinerseits von der Zeit abhängt. Wenn die Durch- 


158 


“ 


senkrechtem Auftreffen der Stoßmasse m dessen ge- 


Be größten Wert d erreicht hat, hat sie bei worinx von Form 


samte Energie aufgenommen und es gilt, wenn von. Dur 


einer Verformungsarbeit am Stoßkörper abgesehen 
wird. u { 


d Y d d 
% 0) 
AN P(o)ds= m Or ds 
91? 01 (4) 
s=0 = 0 \ 8-0 . 7 
| ER ND, | 
he ana 


Die Scherbeanspruchung der Platte wird ursächlich 
von der Verzögerungskraft P der Stoßmasse bestimmt. 
Da dos Integral in Gl.(4) unabhängig von der Kon- 
struktion der Platte konstant, d.h. dem gegebenen 


ER mi ist, wird die,Kraft P um so kleiner, 


je langsamer die Stoßgeschwindigkeit abgebremst wird. 
Da lee daß die Platte möglichst weich und nach- 
giebig zu konstruieren ist, wenn es darauf ankommt, 


die Stoßkraft einer auftreffenden Masse und damit die 


Scherbeanspruchung des Materials herabzusetzen. Die 
Platte wird aber auch auf Biegung beansprucht, die 
strenge Lösung des Schwingungsvorganges beim sym- 
metrischen Biegestoß gegen eine Kreisplatte wurde vom 
Verfasser gegeben, ebenso eine als sehr gut erwiesene 


. Näherungslösung für den Biegepfeil, dieauch auf Recht- 


eckplatten ausgedehnt wurde !). a 
Zur Ausbreitungsgeschwindigkeit der Biegungswel- 
len siehe W. Flügge?). Es soll untersucht werden, 
ob es vorteilhaft ist eine starke Platte in mehrere 
schwächere aufzuteilen, sie also weicher zu machen. 


I. Vergleich verschieden steifer 
Plattenkonstruktionen. 


Um das Wesentliche klar herauszuschälen, werden in 
den folgenden Formeln alle Werte bis auf die Platten- 
masse M, Stoßmasse m, Stoßgeschwindigkeit v,, dem 
statischen bzw. dynamischen Biegepfeil f bzw. d, die 
Zeit T vom Stoßbeginn bis zum Erreichen dieses Biege- 


pfeiles und die Deckenstärke h, in Beiwerte ce zusam- 


mengefaßt. 
Bei statischer Beansprychung ist: 


ea e=@hıf . . ...(da). 
Man sieht sofort, daß bei gleicher Last P eine Auf- 
teilung der Vollplatte, auf die sich der IndexI bezieht, 
in n Platten gleicher Stärke hı/n auf die Last anteilig 
verteilt ist, zu n-facher Biegebeanspruchung führt: 


BERRY N 5 
On==Ca (kım)3 n.07.. (0b), 


BD 
01 = (3 m A 
Wenn aber nicht die statische Last ?, sondern die 
geforderte Arbeitsaufnahme ertragen werden soll, sind 

nach G1.(6) beide Konstruktionen gleichwertig. 

0% Ohr : 
A1=2,.,M> lo Fi = AT% (6). 
Eindeutig zugunsten der weichen Plattenkonstruk- 
tion verschiebt sich das Bild bei dynamischer 
Scherbeanspruchung. Der Stoß wird langsamer ab- 


z ov 
gefangen und damit die der Druckkraft P = m Fr 


proportionale Scherkraft geringer. 
Diese Aussage wird erhärtet durch die zitierte Ar- 
beit 1), wonach als dynamischer Biegepfeil folgt: 


2, m 
m RT 


!) G.Sonntag: Der symmetrische Biegestoß gegen eine Krei- 
selplatte . Forsch, a. d. Gebiet. d. Ing.-Wesens Bd. 14 (1943) 8.137. 
?) W. Flügge: Die Ausbreitung von Biegungswellenin Stäben. 


Z. angew. Math. Mech, Nr. 6 (1942) 8.312. 


der Bi Bnnsce mn ur PLA einander 
a re een 
geschichteten Platte zu vergleichen, Je 


zunächst stillschweigend ein reibungsfreies Aufein- 
anderliegen der Platten v werden muß. Bei 


dem durehzuführenden Vergleich kann daher die Biege 
beanspruchung o proportional d- h ee 
und es folgt t 2, ° 


. . m } 
o=4%\ yır“um) N - 8). 


also unab ig von der Art der Aufteilung der Platte 
in Sehichten Kind die Plattenschichten nicht gekop- 
pelt, sondern soweit distanziert, daß sie nur nachein- 
ander zerstört Tan nach Gl. (8) 
die Beanspruc r Ei in geringerem 
Maße rg ihrer Schwächung gegenüber der 
Vollplatte entspräche, da die Plattemasse M unter der 
Wurzel steht. Die Folgerungen daraus können aber zu 
Trugschlüssen führen, da ein chungsvergleich 
bei verschiedenem Verhältnis m/M zunächst eine um- 


fangreichere Auswertung der zitierten Arbeit *) er- 


fordert. 

Die Scherbeanspruchung r ist bei gegebener Stoß- 
masse und Geschwindigkeit rtional der Stoß- 
verzögerung und der reziproken ‚Plattenstärke. 


Be 7 'ı SR 2 = (9). 


Die maximale Stoßverzögerung b setzen wir propor- 
tional v,/7T. Zur Berechnung von 7 stehen die Formeln 
der Grenzfälle 


m> M a ud m<M T-2 (10) 


zur Verfügung, woraus bei gekoppelten Plattenschich- 
ten folgt: 


m>M n=/mh, ar 
r (11). 


-. 


Die Überlegenheit der weicheren Platte liegt klar zu- 
tage, praktische Bedeutung hat nur der Fallm< M, 
bei dem die Scherbeanspruchung unabhängig von der 
Plattenstärke ist. 

Man könnte einwenden, daß für einen schnellen Stoß 
obige Formeln für Schwingungszeit, Biegepfeil und 
Stoßverzögerung nicht gelten, da der Stoß zur Wirkung 
kommt, bevor die Biegewelle sich bis zum Plattenrand 
ausgebreitet hat. Dem ist zu entgegnen, daß diese 
Größen und besonders die Biegespannungen, die in 
starkem Maße von den Oberschwingungen mitbestimmt 
werden, sich freilich nicht ermitteln lassen; der hier 
beabsichtigte Vergleich jedoch auch während der Aus- 
breitungszeit der Wellen angenähert zutreffen wird, da 
die Laufgeschwindigkeiten der Momenten sowie der 
Querkraftwellen (siehe ’)) unabhängig von der Platten- 
stärke A sind. Deshalb wird angenommen, daß die für 
den Vergleich erforderliche Ähnlichkeit gewahrt bleibt. 


m<M 


€ 
U=-0mM, 
vs 


I. Zusammenfassung. 


Die Vergleichsuntersuchung hat eindeutig ergeben, 
daß eine starke Platte dynamischer Beanspruchung 
weniger Widerstand entgegensetzt, als wenn sie in meh- 
rere schwächere Platten aufgeteilt wird. Diese Tat- 
sache ist meines Erachtens für den Beschuß von Stahl- 
platten schon allgemein bekannt. Wenn man die Zer- 
störung der obersten Platten zulassen will, ist es zweck- 
mäßig, die Platten soweit zu distanzieren, daß genü- 
gend Stoßenergie durch plastische Verformung ver- 
nichtet werden kann. Eine Kopplung der Platten durch 


- 


op - .: 4 
AN ni nt a irn une ey mn m rn nn . 


len der Zwischenräume 2. B. mit Sand erscheint Die Bere Platte sollte öglich ae PR. 
en 2.B.! nd er Markt re Platte sollte möglichst zäh gebaut sein,um 
a onrenenger Biegebeanspruchung durch große  beiihrer Zerstörung einen Be Teil ee A 
Massen kleiner Geschwindigkeiten als zweckmäßig, hin- Arbeitsfähigkeit der auftreffenden Masse zu ver- 


sichtlich der Scherbeanspruchung durch schnelle Mas- schlucken. . 

sen 18 ._ P Be A Fi ” 3 k R n x fr al 
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logarithmisches und nu 
Rechnen. 2. Aufl. 36 $. Bonn 1948. 
lers Verlag. Preis brosch. 2,40 DM. 
. H. Gravelius, Vierstellige Tafeln für 
logarithmisches und numerisches 
Rechnen. Ausgabe A mit Sechsteilung des Grades. 


erisches 
erd. Dümm- 


RER _ BUCHBESPRECHUNGEN 
Prof. Dr. J. Peters, DreistelligeTafelnfür 
 dreisprachigen von Tollhausen ist, nur zwei Sprachen 


aus dem zuerst 1877 im Tauchnitz-Verlag erschienenen R x 


und enthält nur etwa die Hälfte der ursprünglich vor- 
gesehenen Stichworte. Es bringt daher nur die wichtig- 


sten Fachausdrücke der industriellen Technik, wobei 
Hilfswissenschaften, wie Mathematik, Physik und 


Chemie, allerdings nur in ganz beschränktem Umfang, 


Dritte Aufl. 32 8. Bonn 1948. Ferd. Dümmlers Verlag. und das Bauwesen eingeschlossen sind. Weiterer Raum Ba 
Preis brosch. 2,50 DM. x Ä ‚wird dadurch gespart, daß Hauptwörter, die leicht aus Be 
‚. Bei den Tafeln für das logarithmische Rechnen sind Zeit- oder Eigenschaftswörtern abgeleitet werden Be 
in der dreistelligen auch die Antilogarithmen, und die können, nicht aufgeführt sind, daß fast vollständig auf » 
Additions- und Subtraktionslogarithmen gegeben. In tormal gebildete Wortzusammensetzungen sowie auf au 

B der vierstelligen sind sämtliche Tafeln in senkrechten Zusammensetzungen mit Worten gleicher Bedeutung % 

E Spalten angeordnet, in Abweichung z.B. von den verzichtet ist. So ist ein sehr brauchbares Wörterbuch Ki 

3 üblichen Tafeln der Logarithmen. In ihr sind die nicht zu großen Umfanges entstanden. 

= ge der vier Winkelfunktionen von 0° bis 9 Dresden. Willers. 3 

3 miteinem Argumentschritt von 1’, die übrigen miteinen A $ br Ar 

solchen on or gegeben, während die Sdheierälligen Dr. Walter Lietzmann (Prof. a. d. Univ. Göttingen): En 

Sonderlingeim Reich der Zahlen. 1758. a 


3 Be RE To ndche m Er ee en t Sonderbarkeiten haben, und es gibt Zahlen und Zahlen- Re 

E era 1L Tafel bri : “F gegeben. 10 zusammenstellungen, die ihrerseits zu den Sonderlingen = 

2 vierstellige Tafel bringt außerdem Formeln und Er- „anlen; auch sie sind häufiger als man gewöhnlich Fr 

, Iäuterungen zum Gebrauch der Tafel, Konstanten, eine meint.” .Dies kleine Büchlein will diese Menschen und ; 

; ER sche Tnterpolationstabelle usw. Ne . diese Zahlen zusammenbringen“. Beginnend mit In- 

5 resden. Willers. dividuen, die nach seiner Meinung Sonderlinge sind, g 
\8 


Tafeln die Logarithmen aller sechs Winkelfunktionen 
. von 0,1° zu 0,1° aufzeichnen. In den Teilen, die dem 


“numerischen Rechnen dienen sollen, sind in beiden 


Tafeln, neben den üblichen für i/n, n? und in der drei- 


Dr.-Ing. Richard Ernst, Deutsch - Engli- 
sches, Englisch-Deutsches Wörter- 
buch in 2Bänden. Band 1: Deutsch-Englisch 605 8. 
Hamburg 1948. Edition Tauchnitz. Preis geb. 16,50 DM. 


Technische Wörterbücher veralten bei der schnellen 


Entwicklung der Technik bald. Daher ist das Erschei- 


nen eines auch die neuere Entwicklung der Technik 
berücksichtigenden Wörterbuches sehr zu begrüßen. 
Ein die ganzen Ingenieurwissenschaften umfassendes 
Wörterbuch würde sehr umfangreich sein. Um das zu 
vermeiden, umfaßt das vorliegende, das ein Auszug 


m. 5Abb. Bonn 1948. Ferd. Dümmlers Verlag. Preis 


brosch. 6,—DM. 
Verf. schreibt im Vorwort: ‚‚Es gibt Menschen, weit 
mehr als man gemeinhin annimmt, die Freude an 


geht Verf. zu Rechenoperationen weiter, in denen solche 
Sonderlinge auftreten. Der Bereich der Zahlen wird 
dabei von den natürlichen bis zu denreellen ausgedehnt, 
die Operationen von den elementaren _ endlichen 
rationalen zu den unendlichen, die auf unendliche _ 
Dezimalbrüche, Kettenbrüche, Reihen und Produkte 
führen. i 

Das Buch bedeutet eine erfreuliche Bereicherung der 
Bücher, die der Unterhaltungsmathematik gewidmet 
sind, und mancher wird an den einzelnen Abschnitten 
seine Freude haben. 


Dresden. \ Willers. 


Die besprochenen und angezeigten Bücher sind durch den Buchhandel zu beziehen 


EINGEGANGENE BÜCHER 


Bei der Schriftleitung sind folgende Bücher eingegangen (ausführliche Besprechung bleibt vorbehalten) 


Karl Reicheneder, Fehlertheorieund Aus- 
gleichung von Rautenketten in der 
Nadirtriangulation. (Veröffentlichungen des 
Geodätischen Instituts in Potsdam.) 98 S. Akademie- 
Verlag Berlin 1949. Preis brosch. 11,— DM. 


Dr. Heinrich Blasius (Oberstudienrat an der In- 
genieurschule zuHamburg), Wärmelehre. Phy- 
sikalische Grundlagen vom techni- 
schen Standpunkt. 294 S. mit 126 Abb., 
67 Aufgaben und 15 Tabellen. Hamburg 1949. Boysen 
& Maasch Verlag. Preis brosch. 11,50 DM. 


Dr. Heinrich Blasius (Oberstudienrat an.der In- 
genieurschule zu Hamburg), Mechanik. Physi- 
kalischeGrundlagenvomtechnischen 
Standpunkt.. Zweiter Teil: Elastizität 
und Festigkeit. VIII/230 8. m. 194 Abb. und 
84 Aufgaben. Hamburg 1949, Boysen & Maasch Verlag. 
Preis brosch. 9,50 DM. 


Dr. Walter Seheidt (o. ö. Prof. für Anthropologie an 
der Univ. Hamburg), Anthropologie. XI.Zäh- 
lenund Rechnen. 438. m. 34 Abb. Hamburg 
1949. Richard Hermes Verlag. Preis brosch. 2,60 DM. 


Th. Pöschl (o. Prof. a.d. Techn. Hochschule Karls- 
ruhe). Einführung in die Analytisiche 
Mechanik. (Reihe: Mathematische Methoden in 
Naturwissenschaft und Technik, herausgeg. von Dr.» 
Th. Pöschl). ‘166 8. mit 37 Abb. Karlsruhe 1949, 
G. Braun-Verlag. Preis: 'Halbl. gebd. 10,— DM. 


Helmut Hasse, Invariante Kennzeich- 
nung relativ-Abelscher Zahlkörper 
mit vorgegebener Galoisgruppeüber 
einem Teilkörper des Grundkörpers. 
(Abhandlung der Deutschen Akademie der Wissen- 
schaften zu Berlin). 568. Berlin 1949, Akademie- 
Verlag. Preis brosch. 6,— DM. | 


Re - 
. Die Taylor - Entwicklung von YI+z ergibt: vom Verfasser 


Karlsruhe: Der ehem. o. Prof. der Mathematik 
an der Universität Straßburg Dr. Karl Stru- 
becker wurde zum o. Prof. der Mathematik, der 
apl. Prof. an der Universität Göttingen, Dr. Hans 

ittich, zumpl.a. o. Prof. für Mathematik an der 
Techn. Hochschule Karlsruhe ernannt. 


Heidelberg: Am 4. April verschied im Alter 
von 60 Jahren beim Besuch einer Mathematikertagung 
in Oberwolfach der o. Prof. der Mathematik an der i 
Universität Heidelberg, Dr. W. Threlfall. 
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ZUSCHRIFTEN AN DEN HERAUS 


Zu F.K. Rubbert: „Zur Radizierung mit der in etwas anderer 


ei 


Ei 


u 


Rechenmaschine. Z. angew. Math. Mech. 28 ne SE 

(1948) 190—191. Yi+s:-l+ ou 
In Z. angew. Math. Mech.28 Heft 6 gibt Herr F. K. IE EB 

Rubbertin(l) eine Formel, die sich schon bei Vega, 

„‚Logarithmische, trigonometrische uud andere zum 


Gebrauch der Mathematik eingerichtete Tafeln und j 
Formeln“, Wien 1783 8..387 findet. Die in (2) ange- Er gewinnt sie aus einer 


- gebene Fehlergrenze gibt nur die ungefähre Größen- Bild tz für 


ordnung des Fehlers an, 


yi ae eu 
| | en z Ki Er Baalaar Te Beider prakt: r 
LT DE ee | fachere Form des Restgliedes, da man ja bereits von 
+ 6n3 Rh ahe ‚einem relativ guten Näherungswert (vom Rechen 
„Date +nt@n—N@—del  Zano Sellensahl bestimmen mil Re 
bmPpn+Mm—)x] „Feet wirklich enges ch ee RS 
—1 1 2n— I) m—1 urzelformel in Vergessenheit ten konnte, so ut 2rn 
ee Er man sioin keiner Formelsammlungundinkeinem Lehr: 
buch der numerischen Methoden findet, Undin dm 
Seele Arithmetikbuch von Stolz- Gmeiner steht sie 
m ” nur in einem ispi 
o ischen fand der Verfasser — mit elliptischen 
für a ir RE Integralen beschäftigt — eine weitere Literaturstelle, 
3(3n—1) a ee Ch. W. Beste 
o Ä k. ield, OnaNew Method of Rewrepgihe er 
DERRENACN Me to Rllipticand Ultra-elliptie Functions 


Verf. beschreibt ein Verfahren, das mit jedem Schritt Philosophical Transaetions 1862, p. 417-435. H wi 


etwa dreimal so viele Stellen der Wurzelwerte‘liefert wird auch die allgemeine Formel für kubische A i. 
als der vorangehende, Ich darf dazu vielleicht: bemer- mation angeführt: 


ken, daß dieses Verfahren an sich nicht neu, sondern „5 
schon recht alt ist, Es steht bei J. H. Lambert, 
Taytrige zum le a are Lo, En (n2-1)am +(4n?+2)arn R-+(n?-1) R? s 
wendungen, Zweyter Theil, Erster Abschnitt, 8. u re ER TEE TS TEE AT ac HE ET 
(1770) und wird als etnberisikied Verfahren‘‘ gebracht n* SER. SARA te uw. Ian R+n-1)(n-2)R* 
von O. Stolz und J. A. Gmeiner, Theoretische WO a einen ersten Näherungswert bedeutet. Zur Ab- 
Arithmetik IL (2. Aufl. 1915) 8.125, Fußnote 1 und leitung der Nüherungsformel höherer Ordnung wendet 
S. 217, mitsamt der auch von Herrn Rubbert an- Merrifie Id erzeugende Funktionen ann Lam- 
gegebenen Fehlerabschätzung. Ich selbst habe es wie- bert wird nicht genannt, dagegen der bekannte eng- 
derholt in meinen Vorlesungen vorgetragen und freue lische Mathematiker Sylvester, der für n=2 die 
mich, daß es Herr Rubber t nun wieder einem grö- höheren Approximationsformeln gefunden und im Phi- 
Beren Kreis in Erinnerung gebracht hat, von seinen losophical Magazine, vol. 20 (1860) veröffentlicht hat, 


Vorgängern allerdings wohl nichts wissend, Bezeichnet-man die i-te Näherung von / R mit N;, so 
Tübingen, Max Müller. gilt nach Merrifield die allgemeine Formel 
Erwiderung. a ae Ei 
Die Zuschrift von Herrn Professor Max Müller (a-+ Y Rli—(a—YR)! 


hat mich sehr erfreut. Also Lambert war dererste, die stets rational ist. Wie der Titel der genannten Ar- 
der diese so nützliche Wurzelformel fand; ich habe beit: schon sagt; benutzt: ihr Verfasser diese Wurzel- 
gleich in den „‚Beyträgen“‘ nachgelesen, um seine Ab- formeln zur genäherten Quadratur. \ 
leitung kennen zu lernen. Lam berts Formel lautet Göttingen, ? F.K.Rubbert. 
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